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RESEARCH ARTICLE 


Unravelling species boundaries in the Aspergillus viridinutans 
complex (section Fumigati): opportunistic human and animal 
pathogens capable of interspecific hybridization 


V. Hubka'?3*, V. Barrs**, Z. Dudova'**, F. Sklenář, A. Kubatova', T. Matsuzawa, 
T. Yaguchi®, Y. Horie®, A. Novakova?, J.C. Frisvad’, J.J. Talbot*, M. Kolařík? 


Abstract Although Aspergillus fumigatus is the major agent of invasive aspergillosis, an increasing number of 
infections are caused by its cryptic species, especially A. /entulus and the A. viridinutans species complex (AVSC). 
Their identification is clinically relevant because of antifungal drug resistance and refractory infections. Species 
boundaries in the AVSC are unresolved since most species have uniform morphology and produce interspecific 
hybrids in vitro. Clinical and environmental strains from six continents (n = 110) were characterized by DNA se- 
quencing of four to six loci. Biological compatibilities were tested within and between major phylogenetic clades, and 
ascospore morphology was characterised. Species delimitation methods based on the multispecies coalescent model 
(MSC) supported recognition of ten species including one new species. Four species are confirmed opportunistic 
pathogens; A. udagawae followed by A. felis and A. pseudoviridinutans are known from opportunistic human infec- 
tions, while A. felis followed by A. udagawae and A. wyomingensis are agents of feline sino-orbital aspergillosis. 
Recently described human-pathogenic species A. parafelis and A. pseudofelis are synonymized with A. felis and 
an epitype is designated for A. udagawae. Intraspecific mating assay showed that only a few of the heterothallic 
species can readily generate sexual morphs in vitro. Interspecific mating assays revealed that five different species 
combinations were biologically compatible. Hybrid ascospores had atypical surface ornamentation and significantly 
different dimensions compared to parental species. This suggests that species limits in the AVSC are maintained by 
both pre- and post-zygotic barriers and these species display a great potential for rapid adaptation and modulation 
of virulence. This study highlights that a sufficient number of strains representing genetic diversity within a species 
is essential for meaningful species boundaries delimitation in cryptic species complexes. MSC-based delimitation 
methods are robust and suitable tools for evaluation of boundaries between these species. 
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INTRODUCTION hazardous mycotoxins), biotechnology and pharmacology 
(production of bioactive substances, heterologous proteins) (Pitt 
& Hocking 2009, Meyer et al. 2011, Frisvad & Larsen 2015b, 


Sugui et al. 2015, Gautier et al. 2016). 


Aspergillus sect. Fumigati includes approximately 60 species 
occurring predominantly in soil (Hubka et al. 2017). Many are 
of considerable medical importance as they cause human and 
animal infections (Balajee et al. 2005b, 2009, Katz et al. 2005, 
Yaguchi et al. 2007, Hubka et al. 2012, Talbot & Barrs 2018). As- 
pergillus fumigatus is usually reported as both the most common 
member of the section in soil worldwide and the most common 
cause of aspergillosis (Klich 2002, Domsch et al. 2007, Mayr & 
Lass-Fl6orl 2011). A series of recent studies highlighted the high 
prevalence (11—19 %) of so-called cryptic Aspergillus species in 
clinical samples (Balajee et al. 2009, Alastruey-Izquierdo et al. 
2013, Negri et al. 2014, Sabino et al. 2014). Their identification 
is clinically relevant since many demonstrate drug resistance 
to commonly used antifungals, thus their recognition influences 
therapeutic management. Reliable identification of clinical iso- 
lates to the species level and susceptibility testing by reference 
methods is thus warranted (Lyskova et al. 2018). Many of these 


Aspergillus is a speciose genus with almost 400 species classi- 
fied into six subgenera and approximately 25 sections (Samson 
et al. 2014, Jurjević et al. 2015, Hubka et al. 2016a, 2017, Chen 
et al. 2016a, b, 2017, Kocsubé et al. 2016, Sklenář et al. 2017, 
Tanney et al. 2017). The species are widely distributed in nature 
and have a significant economic impact in human and animal 
health (causative agents of aspergillosis; allergies and respira- 
tory problems associated with presence of fungi in the indoor 
environment), the food industry (source of enzymes and organic 
acids for fermentation, food and feed spoilage, production of 


Department of Botany, Faculty of Science, Charles University, Benatska 
2, 128 01 Prague 2, Czech Republic. 

Laboratory of Fungal Genetics and Metabolism, Institute of Microbiology 
of the CAS, v.v.i, Vídeñská 1083, 142 20 Prague 4, Czech Republic. 
First Faculty of Medicine, Charles University, Kateřinská 32, 121 08 Prague 
2, Czech Republic. 

Sydney School of Veterinary Science, Faculty of Science, and Marie Bashir 
Institute of Infectious Diseases & Biosecurity, University of Sydney, Camper- 
down, NSW, Australia. 

University of Nagasaki, 1-1-1 Manabino, Nagayo-cho, Nishi-Sonogi-gun, 


N 


° 


> 


a 


Nagasaki 851-2195, Japan. 

Medical Mycology Research Center, Chiba University, 1-8-1, Inohana, 
Chuo-ku, Chiba 260-8673, Japan. 

Department of Biotechnology and Biomedicine, Technical University of 
Denmark, Kongens Lyngby, Denmark. 

* corresponding author e-mail: hubka@biomed.cas.cz. 

# These co-authors contributed equally to this work. 


D 


s 


less common pathogens belong to sect. Fumigati and the high- 
est numbers of infections are attributed to A. lentulus, A. thermo- 
mutatus (syn. Neosartorya pseudofischeri) and species from 
A. viridinutans species complex (AVSC) (Balajee et al. 2005a, 
2006, Sugui et al. 2010, 2014, Barrs et al. 2013, Talbot & Barrs 
2018). 
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Homothallism is a predominant reproductive mode in sect. 
Fumigati and many species readily produce ascomata (neosar- 
torya-morph) in culture, while others are heterothallic or have 
an unknown sexual morph (Hubka et al. 2017). Homothallic 
species are infrequently pathogenic, although A. thermomutatus 
is a notable exception. The majority of clinically relevant species 
belong to the A. fumigatus clade (Balajee et al. 2005b, 2009, 
Yaguchi et al. 2007, Alcazar-Fuoli et al. 2008) or the AVSC (Su- 
gui et al. 2010, 2014, Barrs et al. 2013, Novakova et al. 2014) 
and are heterothallic. A cryptic sexual cycle of several of these 
opportunistic pathogens, including A. fumigatus (O'Gorman et 
al. 2009), A. /entulus (Swilaiman et al. 2013) and A. felis (Barrs 
etal. 2013), was discovered recently by crossing opposite mat- 
ing type isolates in vitro. 


MAT 1-2-1 (LT796765) 


MAT1-2-1 
MAT1-1-1 (LT796762) 


MAT1-1-1 
MAT 1-2-1 (HF937388) 


MAT 1-2-1 
MAT1-1-1 (LT796763) 
MAT1-2-1 (LT796764) 


MAT1-1-1 
MAT1-1-1 
MAT1-2-1 
MAT1-1-1 
MAT1-2-1 
MAT 1-2-1 


MAT locus*® 
Filamentous fungus collection 


working collection of the Applied and Industrial Mycology 


Collection at the Medical Mycology Research Centre, Chiba Uni- 


Japan Collection of Microorganisms, Tsukuba, Japan; KACC 


Korean Agricultural Culture Collection, 


Agricultural Research Service Culture Collection, Peoria, Illinois, USA. 


Molecular methods are routinely used for identification of spe- 
cies from sect. Fumigati due to overlapping morphological fea- 
tures of their asexual morph. In contrast, the morphology of the 
sexual morph, especially of ascospores, is amongst the most 
informative of phenotypic characteristics in sect. Fumigati. The 
taxonomy of AVSC has developed rapidly since eight of the 
currently 11 recognized species were described in the last four 
years (Barrs et al. 2013, Eamvijarn et al. 2013, Nováková et al. 
2014, Sugui et al. 2014, Matsuzawa et al. 2015, Talbot et al. 
2017). The species boundaries delimitation was usually based 
on comparison of single-gene phylogenies and principles of 
genealogical concordance. In addition, some studies supported 
the species concept by results of in vitro mating experiments 
between opposite mating type strains. With the increasing 
number of species, available isolates and new mating experi- 
ment data, the species boundaries in AVSC became unclear as 
pointed out by Talbot et al. (2017) who used the designation ‘A. 
felis clade’ for A. felis and related species. Importantly, Sugui et 
al. (2014) and Talbot et al. (2017) identified that interpretation 
of in vitro mating assays in sect. Fumigati may be problematic 
because different phylogenetic species in the AVSC were able 
to produce fertile ascomata when crossed between themselves. 
Some even mated successfully with A. fumigatus s.str. 


Culture Collection of Fungi, Prague, Czech Republic; CM 
Czech Republic; DTO 


ëjovice, 


culture collection of the DTU Systems Biology, Lyngby, Denmark; IFM 


Collection of Microscopic Fungi, Institute of Soil Biology, Academy of Sciences of the Czech Republic, Ceské Bud 


CABI's collection of fungi and bacteria, Egham, UK; JCM 
Biological Resource Center, National Institute of Technology and Evaluation, Chiba, Japan; NRRL 


male neutered; ND = not determined. 


Here we presenta critical re-evaluation of species boundaries in 
the AVSC. We examined a large set of clinical and environmen- 
tal strains collected worldwide. We did not use classical phylo- 
genetic methods or genealogical concordance phylogenetic 
species recognition rules (GCPSR) for species delimitation due 
to their unsatisfactory results in previous AVSC studies. Such 
methods, based predominantly on analysis of concatenated 
DNA sequence data or comparison of single-gene phylogenies 
are frequently prone to species over-delimitation or are affected 
by subjective judgements of species boundaries. Instead, we 
used recently introduced delimitation techniques based on coales- 
cent theory and the multispecies coalescent model (MSC) (Flot 
2015). We followed the approach recommended by Carstens 
et al. (2013) that combines species delimitation, species tree 
estimation and species validation steps. Although these meth- 
ods have already been applied to other groups of organisms 
such as animals and plants their use in fungi is scarce (Stewart 
et al. 2014, Singh et al. 2015, Liu et al. 2016, Sklenář et al. 
2017, Hubka et al. 2018). Here, the results of MSC methods 
were taken as a basic hypothesis for species delimitation and 
then further verified by analysis of intra- and interspecific bio- 
logical compatibilities, as well as ascospore dimensions and 
ornamentation. 


Australia, Melbourne, retrobulbar mass in a 1.5-year-old BSH cat, MN, 2010 


USA, Wyoming, Glenrock, soil, mine waste dump, 2010 


Russia, Moscow, soil, < 1968 


USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
China, Urumqi, soil, 2008 


USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
USA, Wyoming, Glenrock, soil, mine waste dump, 2010 
USA, Wyoming, Glenrock, soil, mine waste dump, 2010 


Locality, substrate, year of isolation? 
USA, Wyoming, Glenrock, soil, 2010 


USA, human, clinical material 


CBS culture collection housed at the Westerdijk Institute, Utrecht, The Netherlands; CCF 


Food Fungal Culture Collection, North Ryde, Australia; IBT 


domestic shorthair; FN = female neutered (desexed); MN 


Natural History Museum & Institute, Chiba, Japan; CBS 


CBS 135455 (F29) 


IFM 60854 (F5) 
IFM 60856 (F13) 
IFM 60855 (F31) 


Kasetsart University Fungal Collection, Bangkok, Thailand; NBRC (IFO) 


domestic longhair; DSH 


Belgian Coordinated Collections of Micro-organisms (BCCM/IHEM), Brussels, Belgium; IMI 


CCF 5640 (Yogurt R.) 


CMF ISB 2485 (F12 
CMF ISB 2317 (F10 


CMF ISB 1974 
CMF ISB 2487 (F28 


CMF ISB 1976 
CMF ISB 2495 (F53 


CMF ISB 2491 (F60 


CCF 4383 
IFM 59681 = CCF 4563 
DTO 155-G2 


outgroup 


CMF ISB 2162 


MATERIAL AND METHODS 


(cont.) 
British shorthair; DLH 


Fungal strains 


A. lentulus NRRL 35552! = CBS 117885" = IBT 27201! = KACC 41940" 


CCF 4411 = CMF ISB 1977 


CCF 4170 
CCF 4412 (F9) 
CCF 4413 
CCF 4414 
CCF 4415 
CCF 4416 
CCF 4418 
CCF 4419 
CCF 4420 
IMI 133982 


A total of 110 isolates were examined including new isolates 
and isolates obtained from previously published studies (Katz et 


department housed at the Westerdijk Institute, Utrecht, The Netherlands; FRR 


of the Spanish National Center for Microbiology, Madrid, Spain; CMF ISB 
versity, Japan; IHEM 


Wanju, South Korea; KUFC 
2 Original numbers of strains and personal strain designations are given in parentheses. 


3 BSH 


4 When available, sequence number in public database is given in parentheses; in the remaining cases, the MAT idiomorph was confirmed only on the electrophoretogram (specific PCR and length of amplicons). 


Species / Culture collection nos. +? 
1 Culture collection acronyms: CBM-FA 
5 Sequences generated in this study are in bold. 


A. wyomingensis (cont.) 


Table 1 
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al. 2005, Vinh et al. 2009, Coelho et al. 2011, Shigeyasu et al. 
2012, Barrs et al. 2013, 2014, Eamvijarn et al. 2013, Novakova 
et al. 2014, Sugui et al. 2014, Matsuzawa et al. 2015, Talbot et 
al. 2017) and culture collections. The set comprised 38 clinical 
strains and 72 environmental isolates, including 67 from soil, 
four from cave environments and one from plant material. The 
provenance of isolates is detailed in Table 1. Newly isolated 
strains were deposited into the Culture Collection of Fungi at 
the Department of Botany, Charles University, Prague, Czech 
Republic (CCF). Dried herbarium specimens were deposited 
into the herbaria of the Medical Mycology Research Center, 
Chiba University, Japan (IFM) and Mycological Department of 
the National Museum, Prague, Czech Republic (PRM). 


tsr1 
LT796154 LT796757 


mcm7 


Filamentous fungus collection 


working collection of the Applied and Industrial Mycology 


Collection at the Medical Mycology Research Centre, Chiba Uni- 


Japan Collection of Microorganisms, Tsukuba, Japan; KACC 


Korean Agricultural Culture Collection, 


HF937382 
LT796008 
LT796010 
LT796015 
LT796017 
LT796020 
LT796022 
LT796024 
HF937381 
LT796026 
LT796028 
LT796030 
DQ094860 
LT796032 
LT796036 
DQ094873 KJ914712 KJ914746 


act 


HF937378 
LT796009 
LT796011 

LT796016 
LT796018 
LT796019 
LT796021 

LT796023 
HF937377 
LT796025 
LT796027 
LT796029 
LT796014 
LT796035 
EF669756 


Phenotypic studies 


The strains were grown on malt extract agar (MEA), Czapek 
Yeast Autolysate Agar (CYA), Czapek-Dox agar (CZA), yeast 
extract sucrose agar (YES), CYA supplemented with 20 % 
sucrose (CY20S), and creatine sucrose agar (CREA), and 
incubated at 25 °C. Agar media composition was based on 
that described by Samson et al. (2014). Malt extract and yeast 
extract were obtained from Oxoid (Basingstoke, UK) and Fluka 
Chemie GmbH (Switzerland), respectively. Growth at 42, 45 and 
47 °C was tested on MEA plates sealed with Parafilm. Colour 
determination was performed according to the ISCC-NBS (Inter- 
Society Color Council — National Bureau of Standards) Centroid 
Colour Charts (Kelly 1964). 


Micromorphology was observed on MEA. Lactic acid with cotton 
blue was used as a mounting medium. Photographs were taken 
on an Olympus BX-51 microscope (Olympus DP72 camera) 
using Nomarski contrast. Macromorphology of the colonies was 
documented using a stereomicroscope Olympus SZ61 (with 
Olympus Camedia C-5050 Zoom camera) or Canon EOS 500D. 


Scanning electron microscopy (SEM) was performed using a 
JEOL-6380 LV scanning electron microscope (JEOL Ltd. Tokyo, 
Japan) as described by Hubka et al. (2013b). Briefly, pieces 
of colony or mature ascomata were fixed in osmium tetroxide 
vapours for one wk at 5—10 °C and gold coated using a Bal-Tec 
SCD 050 sputter coater. The specimens were observed using 
40 um spot size and 15—25 kV accelerating voltage. 


Agricultural Research Service Culture Collection, Peoria, Illinois, USA. 


EF669895 


HF933397 
HF933394 
HF933392 
HF933389 
HF933390 
HF933391 
HF933393 
HF933395 
HF933396 
HF933398 
HF933399 
HF933400 
LT796013 
LT796034 


GenBank/ENA/DDBJ accession numbers 
CaM RPB2 


Culture Collection of Fungi, Prague, Czech Republic; CM 


HF933359 
HF933354 
HF933356 
HE578077 
HF933352 
HF933360 
HF933353 
HF933357 
HF933358 
HF933355 
HF933361 
HF933362 
LT796012 
HG426056 HG426053 LT796031 
LT796033 
EF669825 


benA 


HG324081 
HE578062 
EF669969 


ITS 


= culture collection of the DTU Systems Biology, Lyngby, Denmark; IFM 


CABI's collection of fungi and bacteria, Egham, UK; JCM 


Kasetsart University Fungal Collection, Bangkok, Thailand; NBRC (IFO) = Biological Resource Center, National Institute of Technology and Evaluation, Chiba, Japan; NRRL 


Molecular studies 


ArchivePure DNA yeast and Gram2+ kit (5 PRIME Inc., Gaithers- 
burg, MD) was used for DNA isolation from 7-d-old cultures 
according to the manufacturer’s instructions as updated by 
Hubka et al. (2015b). The purity and concentration of extracted 
DNA was evaluated by NanoDrop 1000 Spectrophotometer. 
ITS rDNA region was amplified using forward primers ITS1 or 
ITS5 (White et al. 1990) and reverse primers ITS4S (Kretzer 
et al. 1996) or NL4 (O'Donnell 1993); partial B-tubulin gene 
(benA) using forward primers Bt2a (Glass & Donaldson 1995) or 
Ben2f (Hubka & Kolařík 2012) and reverse primer Bt2b (Glass 
& Donaldson 1995); partial calmodulin gene (CaM) using for- 
ward primers CF1M or CF1L and reverse primer CF4 (Peterson 
2008); partial actin gene (act) using primers ACT-512F and 
ACT-783R (Carbone & Kohn 1999); partial RNA polymerase II 
second largest subunit (RPB2) using forward primers fRPB2-5F 
(Liu et al. 1999) or RPB2-F50-CanAre (Jurjević et al. 2015) and 
reverse primer fRPB2-7cR (Liu et al. 1999); partial mcm7 gene 
encoding minichromosome maintenance factor 7 with primers 
Mcm7-709for and Mcm7-1348rev (Schmitt et al. 2009); and 
partial fsr7 gene encoding ribosome biogenesis protein with 
primers Tsr1-1453for and Tsr1-2308rev (Schmitt et al. 2009). 
Terminal primers were used for sequencing. 


CBS culture collection housed at the Westerdijk Institute, Utrecht, The Netherlands; CCF 
Collection of Microscopic Fungi, Institute of Soil Biology, Academy of Sciences of the Czech Republic, České Budějovice, Czech Republic; DTO 


department housed at the Westerdijk Institute, Utrecht, The Netherlands; FRR 


versity, Japan; IHEM 


Food Fungal Culture Collection, North Ride, Australia; IBT 


IFM 60854 
IFM 60856 
CBS 135455 
IFM 60855 


SB 2317 
SB 1974 
SB 2487 
SB 1976 
SB 2162 
SB 2495 
SB 2491 


CMF ISB 2486 
CMF ISB 2485 
CCF 4383 
CCF 5640 


CMF 
CMF 
CMF 
CMF 
CMF 
CMF 
CMF 


Natural History Museum & Institute, Chiba, Japan; CBS 


NRRL 35552" = CBS 117885" = IBT 27201" = KACC 419407 


CCF 4417" = CMF ISB 2494" = CBS 135456" 


CCF 4169 
CCF 4411 = CMF ISB 1977 


CCF 4412 
IFM 59681 = CCF 4563 


CCF 4413 
CCF 4414 
CCF 4415 
CCF 4416 
CCF 4418 
CCF 4419 
CCF 4420 
IMI 133982 
DTO 155-G2 


Culture collection nos." 
CCF 4170 


Belgian Coordinated Collections of Micro-organisms (BCCM/IHEM), Brussels, Belgium; IMI 


of the Spanish National Center for Microbiology, Madrid, Spain; CMF ISB 


Wanju, South Korea; KUFC 


A. lentulus 


The PCR reaction volume of 20 uL contained 1 uL (50 ng) of 
DNA, 0.3 uL of both primers (25 pM/mL), 0.2 uL of MyTaq™ DNA 


1 Culture collection acronyms: CBM-FA 


Table 2 (cont.) 
Species 
A. wyomingensis 
outgroup 
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Polymerase (Bioline, GmbH, Germany) and 4 uL of 5x MyTaq 
PCR buffer. The ITS rDNA, benA and CaM fragments were 
amplified using the following thermal cycle profile: 93 °C/2 min; 
30 cycles of 93 °C/30 s; 55 °C/30 s; 72 °C/60 s; 72 °C/10 min. 
The annealing temperature for amplification of act gene was 
60 °C (30 cycles); and that for fsr7 gene 50 °C (37 cycles). 
Partial RPB2 gene fragments were amplified using the above- 
mentioned cycle or touchdown thermal-cycling: 93 °C/2 min; 
5 cycles of 93 °C/30 s, 65—60 °C/30 s, 72 °C/60 s; 38 cycles of 
93 °C/30 s, 55 °C/30 s, 72 °C/60 s; 72 °C/10 min. The partial 
mcm7 gene was amplified using modified touchdown thermal- 
cycling: 93 °C/2 min; 5 cycles of 93 °C/30 s, 65—60 °C/30 s, 
72 °C/60 s; 38 cycles of 93 °C/30 s, 60 °C/30 s, 72 °C/60 s; 
72 °C/10 min. PCR product purification followed the protocol of 
Réblova et al. (2016). Automated sequencing was performed at 
Macrogen Sequencing Service (Amsterdam, The Netherlands) 
using both terminal primers. Sequences were deposited into 
the ENA (European Nucleotide Archive) database under the 
accession numbers listed in Table 2. 


Phylogenetic analysis 


Sequences were inspected and assembled using Bioedit v. 7.2.5 
(www.mbio.ncsu.edu/BioEdit/bioedit.html). Alignments of the 
benA, CaM, act and RPB2 regions were performed using the 
G-INS-i option implemented in MAFFT v. 7 (Katoh & Standley 
2013). Alignments were trimmed, concatenated and then ana- 
lysed using Maximum likelihood (ML) and Bayesian inference 
(BI) analyses. Suitable partitioning scheme and substitution 
models (Bayesian information criterion) for analyses were 
selected using the greedy algorithm implemented in Parti- 
tionFinder v. 1.1.1 (Lanfear et al. 2017) with settings allowing 
introns, exons and codon positions to be independent partitions. 
Proposed partitioning schemes and substitution models for each 
dataset are listed in Table 3. The alignment characteristics are 
listed in Table 4. 


The ML tree was constructed with IQ-TREE v. 1.4.4 (Nguyen 
et al. 2015) with nodal support determined by non-parametric 
bootstrapping (BS) with 1000 replicates. Bayesian posterior 
probabilities (PP) were calculated using MrBayes v. 3.2.6 (Ron- 
quist et al. 2012). The analyses ran for 10’ generations, two 
parallel runs with four chains each were used, every 1 000th 
tree was retained, and the first 25 % of trees were discarded as 
burn-in. The trees were rooted with Aspergillus clavatus NRRL 1 
and A. lentulus NRRL 35552, respectively. All alignments 
are available from the Dryad Digital Repository (https://doi. 
org/10.5061/dryad.38889). 


Species delimitation and species tree inference 


Several species delimitation methods were applied to elucidate 
the species boundaries within the AVSC. We followed the 
recommendation of Carstens et al. (2013) and compared the 
results of several different methods. The analysis was divided 
into two parts. Four genetic loci were examined in the first 
analysis which comprised all species from the AVSC while six 
genetic loci were examined in the second analysis focused on 
the clade comprising Aspergillus felis, A. pseudofelis, A. para- 
felis and A. pseudoviridinutans (A. aureolus was used as an 
outgroup). The alignment characteristics are listed in Table 4. 


Only unique nucleotide sequences, selected with DAMBE 
v. 6.4.11 (Xia 2017) were used in the analyses. Nucleotide 
substitution models for particular loci were determined using 
jModeltest v. 2.1.7 (Posada 2008) based on Bayesian informa- 
tion criterion (BIC) and were as follows: 1st analysis - K80+G 
(benA), K80+I (CaM), K80+G (act), K80+G (RPB2); 2nd 
analysis - K80+I (benA), K80+G (CaM), K80 (act), K80 (RPB2), 
HKY+I+G (tsr1), K80 (mcm7). 

In the first analysis, only unique sequences of four loci were 
used, i.e., benA, CaM, act and RPB2. The number of isolates of 
A. felis and A. pseudoviridinutans was reduced to two, because 
this clade was examined in detail in the second analysis based 


Table 3 Partition-merging results and best substitution model for each partition according to Bayesian information criterion (BIC) as proposed by PartitionFinder 


v. 1.1.0. for combined dataset of benA, CaM, act and RPB2 genes. 


Dataset Phylogenetic method 


Partitioning scheme (substitution model) 


Section Fumigati (Fig. 1) Maximum likelihood benA + CaM + act introns (TrNef+G); 3rd codon positions of benA (GTR+G); 1st codon positions of 
benA + CaM + act + RPB2 + 2nd codon positions of act + 3rd codon positions of act (TIM+l); 2nd 


codon positions of benA + CaM + RPB2 (HKY); 3rd codon positions of CaM + RPB2 (HKY+G) 
benA + CaM + act introns (K80+G); 3rd codon positions of benA (GTR+G); 1st codon positions of 
benA + CaM + act + RPB2 + 2nd codon positions of act + 3rd codon positions of act (GTR+1); 2nd 
codon positions of benA + CaM + RPB2 (HKY); 3rd codon positions of CaM + RPB2 (HKY+G) 


Bayesian inference 


A. viridinutans clade (Fig. 5) Maximum likelihood benA + CaM + act introns (K80+G); 3rd codon positions of benA + CaM + RPB2 (TrN+G); 1st codon 
positions of benA + CaM + act + RPB2 + 3rd codon positions of act (TrN); 2nd codon positions of 


benA + CaM + act + RPB2 (F81) 


benA + CaM + act introns (K80+G); 3rd codon positions of benA + CaM + RPB2 (HKY+G); 1st codon 
positions of benA + CaM + act + RPB2 + 3rd codon positions of act (HKY); 2nd codon positions of 
benA + CaM + act + RPB2 (F81) 


Bayesian inference 


Table 4 Overview of alignments characteristics used for phylogenetic analyses. 


Alignment characteristic benA CaM act RPB2 mcm7 tsr1 Combined dataset 
Section Fumigati (Fig. 1) 
Length (bp) 534 697 431 999 _ _ 2661 
Variable position 268 322 234 280 _ _ 1104 
Parsimony informative sites 184 226 148 186 _ _ 744 
A. viridinutans complex (Fig. 5) 
Length (bp) 475 697 344 967 _ _ 2483 
Variable position 115 168 102 135 _ _ 520 
Parsimony informative sites 84 114 70 81 _ _ 349 
A. felis clade (Fig. 3) 
Length (bp) 474 681 329 967 623 761 3835 
Variable position 72 73 35 59 38 103 380 


Parsimony informative sites 50 49 18 32 24 58 231 
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on six loci. Three single-locus species delimitation methods, 
i.e., DGMYC (Reid & Carstens 2012), GMYC (Fujisawa & Barra- 
clough 2013) and PTP (Zhang et al. 2013), and one multilocus 
species delimitation method STACEY (Jones 2017) were used 
to find putative species boundaries. The bGMYC and GMYC 
methods require ultrametric trees as an input, while PTP does 
not. Therefore, single locus ultrametric trees were constructed 
using a Bayesian approach in BEAST v. 2.4.5 (Bouckaert et 
al. 2014) with both Yule and coalescent tree models. We also 
looked at possible differences between strict and relaxed clock 
models, but since these parameters had no effect on the number 
of delimited species, only the results with strict clock model are 
presented here. Chain length for each tree was 1 x 10’ genera- 
tions with 25 % burn-in. The highest credibility tree was used for 
the GMYC method and 100 trees randomly sampled throughout 
the analysis were used for the bGMYC method. Both methods 
were performed in R v. 3.3.4 (R Core Team 2015) using bgmyc 
(Reid & Carstens 2012) and splits (SPecies’ Llmits by Threshold 
Statistics) (Fujisawa & Barraclough 2013) packages. The non- 
ultrametric trees for the PTP method were constructed using 
the ML approach in RAXML v. 7.7.1 (Stamatakis et al. 2008) 
and IQ-TREE v. 1.5.3 (Nguyen et al. 2015) with 1000 bootstrap 
replicates. The PTP method was performed on the web server 
http://mptp.h-its.org/ (Kapli et al. 2017) with p-value set to 0.001. 
The multilocus species delimitation was performed in BEAST 
v. 2.4.5 with add-on STACEY v. 1.2.2 (Jones 2017). The chain 
length was set to 5 x 10° generations, priors were set as fol- 
lows: the species tree prior was set to the Yule model, growth 
rate prior was set to lognormal distribution (M = 5, S = 2), clock 
rate priors for all loci were set to lognormal distribution (M = 0, 
S = 1), PopPriorScale prior was set to lognormal distribution 
(M = -7, S = 2) and relativeDeathRate prior was set to beta 
distribution (a = 1, B = 1000). The output was processed with 
SpeciesDelimitationAnalyzer (Jones 2017). 


The species tree was inferred using *BEAST (Heled & Drummond 
2010) implemented in BEAST v. 2.4.5. The isolates were as- 
signed to a putative species according to the results of the 
above-mentioned species delimitation methods. The MCMC 
analysis ran for 1 x 10° generations, 25 % of trees were discard- 
ed as a burn-in. The strict molecular clock was chosen for all 
loci and population function was set as constant. Convergence 
was assessed by examining the likelihood plots in Tracer v. 1.6 
(http: //tree.bio.ed.ac.uk/software/tracer). We also constructed 
the phylogenetic tree based on concatenated alignment of all 
four loci in |Q-TREE v. 1.5.3 with 1000 bootstrap replicates and 
the optimal partitioning scheme determined by PartitionFinder 
v. 2.1.1 (Lanfear et al. 2017). 


The validation of the species hypotheses was performed in 
BP&P v. 3.3 (Bayesian phylogenetics and phylogeography) 
(Yang & Rannala 2010). The isolates were assigned to the 
species based on the results of species delimitation methods 
and the species tree inferred with *BEAST was used as a guide 
tree. Three different combinations of the prior distributions of 
the parameters 8 (ancestral population size) and T, (root age) 
were tested as proposed by Leaché & Fujita (2010), i.e., large 
ancestral population sizes and deep divergence: 8 ~ G (1, 10) 
and r, ~ G (1, 10); small ancestral population sizes and shal- 
low divergences among species: 8 ~ G (2, 2000) and T, ~ G 
(2, 2000); large ancestral populations sizes and shallow di- 
vergences among species: 8 ~ G (1, 10) and r, ~ G (2, 2000). 
The second analysis with six protein-coding loci, i.e., benA, 
CaM, act, RPB2, mcm7 and tsr1, consisted of the same steps 
as described above. Instead of PTP, we used the programme 
mPTP (Kapli et al. 2017) with IQ-TREE and RAxML trees as 
an input. Within the mPTP programme we used the following 
settings: Maximum likelihood species delimitation inference 
(option ML) and a different coalescent rate for each delimited 


species (option multi). R package ggtree (Yu et al. 2017) and 
the programme densitree (Bouckaert 2010) were used for 
visualization of the phylogenetic trees. 


Mating experiments 


The MAT idiomorph was determined using the primer pairs 
alpha1 and alpha2 located in MAT1-1-1 locus (alpha box 
domain), and HMG1 and HMG2 primers located in MAT 1-2-1 
locus (high-mobility-group domain) as described by Sugui et 
al. (2010). The MAT idiomorphs were differentiated based on 
the different lengths of PCR products visualized by gel electro- 
phoresis; absence of opposite MAT idiomorph was also verified 
in all isolates. The identity of PCR products was proved by 
DNA sequencing in several isolates (accession numbers in 
Table 1); product purification and sequencing were performed 
at Macrogen Europe (Amsterdam, The Netherlands) using 
terminal primers. Selected opposite mating type strains were 
paired within and between major phylogenetic clades on MEA 
and oatmeal agar (OA; Difco, La Ponte de Claix, France) plates 
and incubated at 25, 30 and 37 °C in the dark. The plates 
were sealed with Parafilm and examined weekly from the third 
wk of cultivation for two months under a stereomicroscope 
for the production of ascomata. The presence of ascospores 
was determined using light microscopy. Width and height of 
ascospores were recorded at least 35 times for each success- 
ful mating pair. 


Statistical analysis 


Statistical differences in the width and height of the asco- 
spores of particular species and interspecific hybrids were 
tested with one-way ANOVA followed by Tukey’s HSD (honest 
significant difference) post hoc test in R v. 3.3.4 (R Core Team 
2015). R package multcomp (Hothorn et al. 2008) was used 
for the calculation and package ggplot2 (Wickham 2009) for 
visualization of the results. 


Exometabolite analysis 


The extracts were prepared according to Houbraken et al. 
(2012). High-performance liquid chromatography with diode- 
array detection was performed according to Frisvad & Thrane 
(1987, 1993) as updated by Nielsen et al. (Nielsen et al. 2011). 
Fungi were incubated for 1 wk at 25 °C in darkness on CYAand 
yeast extract sucrose (YES) agars for exometabolite analysis. 


RESULTS 


Phylogenetic definition of AVSC 


In the phylogenetic analysis, 76 combined benA, CaM, act and 
RPB2 sequences were assessed for members of sect. Fumigati. 
The analysis was based on the modified alignment previously 
used by Hubka et al. (2017) and enriched by taxa from AVSC. 
In the Bayesian tree shown in Fig. 1, members of sect. Fumi- 
gati are resolved in several monophyletic clades. The analysis 
showed thatAVSC is a phylogenetically well-defined group and 
the clade gained full support. Similarly, some other clades are 
well-supported by both BI and ML analyses including A. spino- 
sus Clade, A. brevipes clade, A. tatenoi clade, A. thermomutatus 
clade and A. fennelliae clade; A. spathulatus forms a single- 
species lineage distantly related to other clades. Other clades 
have moderate or low support and the species represented 
therein may differ based on genetic loci used for phylogenetic 
reconstruction and taxa included in the analysis. Heterothallic 
species are dispersed across sect. Fumigati (Fig. 1) but the 
majority of them cluster in AVSC and A. fumigatus clades. These 
two clades also encompass the highest number of human and 
animal pathogens in sect. Fumigati not only in terms of their 
number but also their clinical relevance. 
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Reproductive strategy 
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Fig. 1 Phylogenetic relationships of the sect. Fumigati members inferred from Bayesian analysis of the combined, 4-gene dataset of B-tubulin (benA), 
calmodulin (CaM), actin (act) and RNA polymerase II second largest subunit (RPB2) genes. Bayesian posterior probabilities (PP) and Maximum likelihood 
bootstrap supports (BS) are appended to nodes; only PP > 95 % and BS > 70 % are shown; lower supports are indicated with a hyphen, whereas asterisks 
indicate full support (1.00 PP or 100 % BS); ex-type strains are designated by a superscript T, species names in quotes are considered synonyms; the bar 
indicates the number of substitutions per site. The tree is rooted with Aspergillus clavatus NRRL 1°. The reproductive mode of each species is designated by 


icons before the species name (see legend). 
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Fig.2 Schematic representation of results of species delimitation methods in Aspergillus viridinutans species complex based on four genetic loci. The results 
of multilocus method (STACEY) are compared to results of single-locus methods (PTP, BGMYC, GMYC). The results of STACEY are shown as tree branches 
with different colours, while the results of single-locus methods are depicted with coloured bars highlighting congruence across methods. The displayed tree is 
derived from IQ-TREE analysis based on a concatenated dataset and is used solely for the comprehensive presentation of the results from different methods. 
The species validation analysis results (BP&P) are appended to nodes and shown in grey bordered boxes; the values represent posterior probabilities calcu- 
lated in three scenarios having different prior distributions of parameters 8 (ancestral population size) and 1, (root age). The top value represents the results 
of analysis with large ancestral population sizes and deep divergence: 8 ~ G (1, 10) and T, ~ G (1, 10); the middle value represents the results of analysis with 
large ancestral populations sizes and shallow divergences among species: 9 ~ G (1, 10) and r, ~ G (2, 2000); and the bottom value small ancestral population 
sizes and shallow divergences among species: 8 ~ G (2, 2000) and T, ~ G (2, 2000). 
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Fig. 3 Schematic representation of results of species delimitation methods in Aspergillus felis clade based on six genetic loci. The results of multilocus 
method (STACEY) are compared to results of single-locus methods (mPTP, bGMYC, GMYC). The results of STACEY are shown as tree branches with different 
colours, while the results of single-locus methods are depicted with coloured bars highlighting congruence across methods. The displayed tree is derived from 
IQ-TREE analysis based on a concatenated dataset and is used solely for the comprehensive presentation of the results from different methods. The species 
validation analysis results (BP&P) are appended to nodes and shown in grey bordered boxes; the values represent posterior probabilities calculated in three 
scenarios having different prior distributions of parameters 8 (ancestral population size) and T, (root age). The top value represents the results of analysis with 
large ancestral population sizes and deep divergence: 8 ~ G (1, 10) and T, ~ G (1, 10); the middle value represents the results of analysis with large ancestral 
populations sizes and shallow divergences among species: 8 ~ G (1, 10) and 1, ~ G (2, 2000); and the bottom value small ancestral population sizes and 
shallow divergences among species: 8 ~ G (2, 2000) and T, ~ G (2, 2000). 
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Species delimitation and validation in AVSC 


In the first analysis, four genetic loci were examined across 
species of AVSC, isolates of A. felis and its close relatives 
were reduced to two individuals, because a separate analysis 
based on six loci was performed for this clade. Eleven tentative 
species were delimited in AVSC using STACEY. The results 
are summarised in Fig. 2, the differences in the colour of the 
tree branches reflect species delimited by the analysis. The 
analysis supported recognition of three putative species in 
A. udagawae lineage, delimitation of A. acrensis (described 
below) from A. aureolus was not supported, other AVSC spe- 
cies were supported by STACEY without differences from their 
current concept. 


The results derived from STACEY were compared to those 
from three single-locus species delimitation methods. The con- 
sensual results from single-locus species delimitation methods 
are generally in agreement with the results of STACEY for the 
majority of species but vary greatly for A. udagawae, A. aureolus 
and A. acrensis lineages (Fig. 2). Recognition of three putative 
species in A. udagawae lineage was supported only based on 
the CaM locus, while based on benA locus, none of these three 
sublineages gained support. Various delimitation schemes were 
proposed by different single-locus species delimitation methods 
in the A. udagawae lineage based on the RPB2 gene (results 
even varied between the analyses based on different input trees 
for the PTP and GMYC methods), while five putative species 
were identically delimited based on the act locus. The methods 
relatively consistently supported delimitation of the A. acrensis 
lineage based on the RPB2 locus and similarly, BG MYC and 
GMYC methods supported this species based on the act locus. 
In contrast, lineages of A. acrensis and A. aureolus were not 
split by any method when analyzing benA and CaM loci. 


The species validation analysis results are appended to nodes 
of the tree in Fig. 2. Areasonable support is defined by posterior 
probabilities 2 0.95 under all three scenarios simulated by dif- 
ferent prior distributions of parameters 6 (ancestral population 
size) and T, (root age). Delimitation of all putative species (those 
delimited by STACEY, A. acrensis and A. aureolus) were sup- 
ported by the posterior probability 0.98 or higher based on the 
analysis in BP&P v. 3.1 (Yang & Rannala 2010) under all three 
scenarios. The only exception was lower support for splitting of 
A. acrensis and A. aureolus; this scenario was supported by the 
posterior probabilities 0.84, 0.88, 1.00, respectively. 


Species delimitation and validation in A. felis clade and 
its relatives 


In the second analysis, six genetic loci were examined across 
isolates of A. felis, A. parafelis, A. pseudofelis and A. pseudo- 
viridinutans. Only two tentative species, A. felis and A. pseudo- 
viridinutans, were delimited in this clade using STACEY. The 
results are shown as branches designated by different colours in 
Fig. 3. The analysis did not support separation of A. pseudofelis 
and A. parafelis from A. felis; A. fumigatus var. sclerotiorum is 
included in the lineage of A. pseudoviridinutans. 


The results of three single-locus species delimitation methods 
were compared to those from STACEY, and the consensual 
results showed a general agreement (Fig. 3). Delimitation of 
A. pseudofelis from A. felis was not supported by any of the 
used methods. Only a negligible number of analyses supported 
delimitation of basal clades in A. felis as tentative species (des- 
ignated as clade 2 and 3 in Fig. 3). But even in these minority 
scenarios, there were no clear consensual delimitation patterns 
that would support delimitation of A. parafelis. Interestingly, 
mPTP analysis based on act, benA, CaM (with RAxML trees as 
an input only), mcm7 and tsr7 loci together with GMYC analysis 
based on benA (only input tree based on coalescent tree model) 


and act (only input tree based on Yule tree model) loci did not 
support delimitation of A. pseudoviridinutans from a robust clade 
of A. felis. An incomplete lineage sorting was observed between 
A. felis and A. pseudoviridinutans (Fig. 3) evidencing that there 
was probably an ancestral gene flow between these lineages. 
Two isolates from A. felis lineage (IFM 59564 and CCF 5610) 
have benA sequences that cluster with A. pseudoviridinutans 
while sequences of the remaining 5 loci placed them in the 
A. felis lineage (single-gene trees not shown). 


The species validation analysis results are appended to nodes 
of the tree in Fig. 3. Delimitation of A. felis and A. pseudo- 
viridinutans gained absolute support in BP&P analysis (Yang 
& Rannala 2010) under all three scenarios simulated by dif- 
ferent prior distributions of parameters 6 (ancestral population 
size) and 1, (root age). Delimitation of three putative species 
within A. felis lineage gained no support (posterior probability 
0.51) under the scenario with small ancestral population sizes 
and shallow divergences among species: 8 ~ G (2, 2000) and 
T, ~ G (2, 2000). 


Species tree 


The species tree topology was inferred with *BEAST (Heled & 
Drummond 2010) and is shown in Fig. 4. It was used as a guide 
tree during species validation using BP&P but it also represents 
the most probable evolutionary relationships between species 
in the AVSC. The analysis confirmed recombination between 
three subclades of A. felis (Fig. 4) which include also recently 
proposed species A. parafelis and A. pseudofelis thus repre- 
senting the synonyms of A. felis. Similarly, the recombination 
between three subclades of A. udagawae rejected the hypoth- 
esis that they could be considered separate species (Fig. 4). 
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Fig.4 Species tree inferred with *BEAST visualized by using DensiTree 
(Bouckaert 2010). All trees created in the analysis (except 25 % burn-in 
phase) are displayed on the left side. Trees with the most common topology 
are highlighted by blue, trees with the second most common topology by 
red, trees with the third most common topology by pale green and all other 
trees by dark green. On the right side, the consensus trees of the three most 
common topologies are displayed. 
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Fig. 6 Schematic depiction of results of intraspecific mating experiments between opposite mating type isolates of heterothallic members of the Aspergil- 
lus viridinutans species complex. Only successful mating experiments are displayed by connecting lines between opposite mating type isolates; remaining 
mating experiments were negative. Isolates marked by asterisk were only crossed with ex-type strains of A. felis (CBS 130245"), A. parafelis (NRRL 62900") 
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The remaining species delimited in previous steps (Fig. 4), 
including A. pseudoviridinutans and A. acrensis (introduced in 
this study), were supported by *BEAST analysis. The species 
tree had identical topology with the trees inferred by ML and 
BI analyses of the concatenated and partitioned dataset (Fig. 
5), and all species supported by *BEAST had 100 % ML boot- 
strap support (ML BS) and 1.00 BI posterior probabilities (BI 
PP). Several deep nodes in the species tree had only limited 
support similarly to ML and BI analyses. Thus clear positions 
of A. wyomingensis and A. siamensis within the clade also 
containing A. udagawae, A. acrensis and A. aureolus remains 
unresolved, while A. acrensis with A. aureolus form a sister 
clade to A. udagawae (this topology gained absolute support 
in all further analyses — see below). Another robust clade con- 
tained sister species A. felis and A. pseudoviridinutans. The 
remaining species, i.e., A. viridinutans, A. frankstonensis and 
A. arcoverdensis, formed a basal clade in the AVSC and their 
positions within the clade are fully resolved (Fig. 4). 


Clustering of isolates by origin and mating-type idiomorph 


In the phylogenetic analysis, 111 combined benA, CaM, act and 
RPB2 sequences were assessed for members of AVSC. All 
species delimited by methods based on the coalescent model 
were fully supported by BI and ML analyses (Fig. 5). 


The A. udagawae lineage included 25 isolates that clustered 
in three main clades. Mating type gene idiomorph MAT 1-1-1 
was detected in 10 isolates while 14 strains had MAT 1-2-1 idio- 
morph (MAT idiomorph was not determined in one strain). The 
majority of North American isolates (10/14) clustered in clade 1 
together with one strain from Australia; clade 2 comprised only 
three strains originating from Asia; isolates from four different 
continents were present in clade 3. There was no apparent 
clustering based on clinical or environmental origin of strains, 
or their MAT idiomorph. All three clinical isolates from Asia had 
an identical haplotype based on four studied protein-coding 
loci (Fig. 5) but one strain had MAT1-2-1 idiomorph in contrast 
to MAT 1-1-1 idiomorph detected in the remaining two strains. 


The A. acrensis lineage included five strains isolated from soil 
(Brazil) or cave sediment (Romania), two of which had MAT1-1-1 
idiomorph and three had MAT 1-2-1 idiomorph. This lineage is 
very closely related to a homothallic species A. aureolus rep- 
resented by four strains in our analysis. The only known clini- 
cal isolate of A. aureolus (IHEM 22515) was isolated from the 
cornea of a patient in Peru. We were unable to source further 
information about this case and thus the clinical relevance of 
this isolate cannot be confirmed. 


The mutual phylogenetic position of homothallic A. siamensis 
and heterothallic A. wyomingensis remains unresolved. Asper- 
gillus siamensis was represented in our analysis by only two 
isolates from soil in Thailand, which were included in the original 
description (Eamvijarn et al. 2013). The A. wyomingensis line- 
age included 15 isolates; 12 of them came from Wyoming (USA) 
and were closely related to each other and to one isolate from 
China, while two isolates from Australia and Europe displayed 
a higher number of unique positions. The ratio of MAT1-1-1 
isolates to MAT1-2-1 isolates was 8:7, and the majority of 
MAT 1-1-1 isolates from the USA (6/7) clustered in a separate 
subclade that was only supported in the BI analysis. 


The A. felis lineage comprised 35 isolates that clustered in three 
main clades. Mating type gene idiomorph MAT1-1-1 was de- 
tected in 12 isolates, while 20 strains had MAT1-2-1 idiomorph 
(MAT idiomorph was not determined in three strains). There was 
no clustering based on geographic origin as all three clades 
included isolates from two to four continents. Clade 3 contained 
only clinical isolates (n = 4). Clinical strains were predominant 
in clade 1 (18: 2) whereas environmental strains dominated 


in clade 2 (7 : 4). The ratio of MAT1-1-1 isolates to MAT1-2-1 
isolates in clade 1 was balanced (10 : 9) but was biased toward 
MAT 1-2-1 idiomorph in clades 2 (1 : 7) and 3 (1 : 3). Eight iso- 
lates of A. pseudoviridinutans, a sister species of A. felis, were 
examined in this study; MAT1-1-1 idiomorph was determined in 
five of them and MAT1-2-1 idiomorph in three of them. There 
was no apparent clustering based on clinical or environmental 
origin of strains, or their MAT idiomorph (Fig. 5). 


A basal clade of AVSC comprises three soil-borne species. 
Whilst A. viridinutans and A. frankstonensis are known only from 
one locality in Australia, 13 A. arcoverdensis strains included 
in the analysis were isolated on three continents, i.e., South 
America, Asia and Australia. Both, A. viridinutans and A. frank- 
stonensis were represented only by one and two isolates, re- 
spectively, included in the original descriptions (McLennan etal. 
1954, Talbot et al. 2017), and only isolates of one mating type 
are known for each of these species. Isolates of both mating 
types were present in A. arcoverdensis (MAT1-1-1 : MAT1-2-1 
ratio, 8 : 5). Ageographical clustering was apparent in A. arco- 
verdensis strains; two strains from China and one strain from 
Australia formed sublineages separate from the Brazilian 
strains (Fig. 5). 


Mating experiments and morphology of spores 


The MAT1-1-1 and MAT 1-2-1 idiomorphs were determined for 
100 of 104 isolates representing heterothallic species in AVSC 
(Table 1). Systematic mating experiments were first performed 
within major phylogenetic clades of the AVSC. Opposite mating 
type strains representing genetic and geographic diversity for 
each heterothallic species were selected for mating experi- 
ments and crossed in all possible combinations if not otherwise 
indicated (Fig. 6). Successful mating was observed in lineages 
of A. felis, A. udagawae and A. wyomingensis. At least some 
individuals representing all three phylogenetic subclades of 
A. felis (Fig. 3, 5) and A. udagawae (Fig. 2, 5) crossed suc- 
cessfully with individuals from the other subclades. The mating 
capacity of individual isolates was unequal. Whilst some isolates 
of a particular species were able to mate with a broad spectrum 
of opposite mating type strains of the same species, others 
produced fertile ascomata with only a limited set of strains 
or did not mate at all. The morphology of ascospores among 
different crosses in these three species was consistent (Fig. 
7). The exception was great variability in the convex surface 
ornamentation of A. wyomingensis ascospores among and as 
well as within pairings of different isolates ranging from almost 
smooth, tuberculate to echinulate (Fig. 7). Although both the 
width and height of ascospores of A. felis, A. udagawae and 
A. wyomingensis overlapped significantly, their dimensions 
were statistically significantly different (Tukey’s HSD test, 
p value < 0.05) (Fig. 6). No fertile cleistothecia were produced 
by crossing opposite mating type isolates of A. pseudoviridi- 
nutans, A. acrensis and A. arcoverdensis. Mating experiments 
were not performed in A. viridinutans and A. frankstonensis due 
to the absence of opposite mating type strains. 


Opposite mating type isolates of each heterothallic species were 
also selected for interspecific mating assays and crossed in all 
possible combinations. Morphological characteristics of AVSC 
ascospores and induced hybrids are summarised in Table 5. 
Only three of 12 selected A. udagawae isolates produced fertile 
ascomata with some isolates of A. felis, A. wyomingensis or 
A. acrensis (Fig. 8). The highest mating capacity was observed 
in the ex-type strain of A. udagawae IFM 46972 that produced 
fertile ascomata when crossed with A. felis (CCF 5609, CCF 
4171 and CCF 5611), A. wyomingensis (CCF 4416 and CCF 
4411) and A. acrensis (IFM 57290). The width and height of 
ascospores of interspecific hybrids between A. udagawae 
and A. acrensis were significantly different (Tukey’s HSD test, 
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Fig. 7 Comparison of morphology of sexual morphs of A. felis, A. udagawae and A. wyomingensis. a. Fertile cleistothecia of A. felis as a result of crossing 
of isolates IFM 60053 x FRR 5680; b. ascospores in light microscopy; c—d. ascospores in scanning electron microscopy: CBS 130245" x CCF 5627 (c), CBS 
130245" x IFM 60053 (d); e. fertile cleistothecia of A. udagawae as a result of crossing of isolates IFM 469721 x IFM 46973; f. ascospores in light microscopy; 
g—h. ascospores in scanning electron microscopy; i. fertile cleistothecia of A. wyomingensis as a result of crossing of isolates CCF 4416 x CCF 44177; j. 
ascospores in light microscopy (CCF 4416 x CCF 4169); k-n. ascospores in scanning electron microscopy: CCF 4416 x CCF 44177 (k-I), CCF 44177 x CCF 
4419 (m-n). — Scale bars: b, f, j = 5 um; c—d, g—h, k-n = 2 um. 
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Table 5 Ascospores characteristics of Aspergillus viridinutans complex species and interspecific hybrids. 


Ascospore body 


Species / interspecific hybrid (x) (mean + standard deviation; um) 


Ornamentation of ascospores 


width height length of crests (um)' surface ornamentation? 
Aspergillus aureolus 4.8 +0.5 4.4 + 0.4 (0.5-)1-1.5 crests present’; CS tuberculate to echinulate (SEM) 
A. felis 4440.5 3.9 + 0.6 0.5-1.5(—2) crests present’; CS tuberculate to echinulate (SEM) 
A. siamensis 4.5 £ 0:5 3.7+40.4 (0.5-)1-1.5 crests present’; CS tuberculate, echinulate to reticulate (SEM) 
A. udagawae 48+0.4 4.2 + 0.4 (0—)0.5(-1) visible crests absent in < 10 % of ascospores (LM); 
CS tuberculate to reticulate (SEM) 
A. wyomingensis 4.2 + 0.4 3.4 + 0.4 visible crests absent in > 50 % of ascospores (LM); 
CS almost smooth, tuberculate, echinulate (SEM) 
A. felis x A. pseudoviridinutans 4.9 + 0.4 4.2 + 0.5 (0—)0.5-1 visible crests absent in 5—20 % of ascospores (LM) depending 
on parental isolates; CS tuberculate to echinulate (SEM) 
A. felis x A. wyomingensis 4.8 + 0.5 4.3 + 0.3 (0-)0.5-1 visible crests absent in ~ 10 % of ascospores (LM); 
CS tuberculate (SEM) 
A. felis x A. udagawae 5.1+0.5 4.5 + 0.5 0-0.5(-1) visible crests absent in ~ 20 % of ascospores (LM); 
CS echinulate, tuberculate to reticulate (SEM) 
A. udagawae x A. wyomingensis 5.0 + 0.4 4.6+0.3 visible crests absent in ~ 15 % of ascospores (LM); 
CS tuberculate (SEM) 
A. udagawae x A. acrensis 5.2 £ 0.5 4.4 + 0.5 visible crests absent in — 50 % of ascospores (LM) in CMF 


ISB 2190 x IFM 57290 and in 100 % of ascospores in 
IFM 46972 x IFM 57290; CS tuberculate to echinulate in CMF 
ISB 2190 x IFM 57290 (SEM) and glabrous in IFM 46972 x 
IFM 57290 (LM) 


‘Values in parentheses are less common (less than 10 % of measurements). 
2 LM = light microscopy; SEM = scanning electron microscopy; CS = convex surface. 
3 Crests may absent in < 1 % of ascospores in some isolates / crosses. 


p value < 0.05) from A. udagawae (Fig. 8). Approximately 
50 % of hybrid ascospores from mating CMF ISB 2190 with 
IFM 57290 lacked visible equatorial crests and if present, they 
were frequently interrupted or stellate (Fig. 9) in contrast to 
A. udagawae (visible crests present in > 90 % of ascospores, 
crests continuous). The ascomata from mating IFM 46972 
with IFM 57290 contained only low numbers of ascospores 
that were globose or subglobose and glabrous (without crests 
and ornamentation on convex surface). This observation sup- 
ported the hypothesis that A. acrensis is a separate species 
despite its close phylogenetic relationships to A. udagawae. 
The ascospore dimensions of hybrids between A. udagawae 
and A. wyomingensis were similar to those of A. udagawae 
and both width and height were significantly different (Tukey’s 
HSD test, p value < 0.05) from A. wyomingensis (Fig. 8). These 
hybrid ascospores had well-defined equatorial crests that were 
most commonly 0.5—1 um broad and similar to those of A. uda- 
gawae (Fig. 9). The hybrids of A. udagawae and A. felis had 
ascospores with similar equatorial crest length and body width 
to A. udagawae but were significantly different from A. felis, 
and their height was significantly different from both A. felis 
and A. udagawae (Fig. 8). The ascomata of hybrids between 
A. udagawae with A. wyomingensis and A. felis, respectively, 
usually contained only low numbers of ascospores. No mating 
or production of infertile ascomata only was observed between 
crosses of A. udagawae and the remaining heterothallic AVSC 
members (Fig. 8). Interestingly, the majority of interspecific hy- 
brids produced approximately 1—10 % of globose or subglobose 
asco-spores with abnormally large dimensions, approximately 
6.5-10.5 um diam (their dimensions were not included for 
calculations of statistical measures in Fig. 8 and 10, and in 
Table 5). These cells had thick walls similar to normal asco- 
spores, but lacked equatorial crests and had a glabrous or echi- 
nulate surface. Their dimensions were intermediate between 
normal ascospores and asci but their walls were dissimilar to 
those of thin-walled asci. These cells were not observed among 
progeny of the intraspecific crosses (intraspecific mating assay) 
and their presence probably indicates a defect in meiosis and 
ascospore development. 


Two MAT1-1-1 isolates of A. pseudoviridinutans selected for 
interspecific mating assays, namely the ex-type strain NRRL 
62904 and strain IFM 59502, were able to mate with a relatively 
high number of MAT1-2-1 isolates of A. felis (Fig. 10). The 
ascospores of these hybrids were statistically significantly dif- 
ferent in their width and height from A. felis. Equatorial crests 
were absent in approximately 5—20 % hybrid ascospores and, 
if present, they were shorter than those of A. felis (Table 5). 
These observations suggest that A. pseudoviridinutans should 
be treated as a separate species as proposed by species delimi- 
tation methods despite the close phylogenetic relationships of 
both species and incomplete lineage sorting detected between 
these two species (Fig. 3). Only one interspecific hybrid was 
induced in our assay between A. wyomingensis CCF 4169 and 
A. felis NRRL 62900. The ascospore body width and height of 
this hybrid was significantly different from both parental species 
(Fig. 10). In contrast to A. wyomingensis, equatorial crests were 
present in the majority of hybrids and they were occasionally 
interrupted and stellate (Fig. 11). Infertile ascomata were ob- 
served in some crosses between A. felis and following species: 
A. acrensis, A. wyomingensis and A. viridinutans. 


Aspergillus aureolus and A. siamensis are the only two homo- 
thallic species in the AVSC and readily produce ascomata ona 
broad spectrum of media and growth temperatures and are eas- 
ily distinguishable from the eight heterothallic AVSC members. 
Most A. aureolus isolates produce distinctive yellow colonies in 
contrast to the whitish colonies of A. siamensis (Fig. 12). The 
ascospores of both species have similar dimensions, convex 
surface ornamentation and equatorial crest length (Table 5, 
Fig. 12) and most closely resemble those of A. felis from among 
heterothallic species. 


The macromorphology of colonies, micromorphology of asexual 
morphs and physiology have only limited discriminatory power 
in AVSC members, as recognized in previous studies (Novakova 
et al. 2014, Matsuzawa et al. 2015). We compared surface 
ornamentation of conidia in all currently recognized species 
using SEM. The ornamentation showed a micro-tuberculate 
pattern and was broadly identical across all species (Fig. 13). 
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Fig.8 Schematic depiction of results of interspecific mating experiments between opposite mating type isolates of A. udagawae and other heterothallic 
members of Aspergillus viridinutans species complex. Only successful mating experiments are displayed by coloured connecting lines between opposite 
mating type isolates (different colours correspond to hybrids between different species); grey dashed lines indicate production of infertile ascomata; remain- 
ing mating experiments were negative. Boxplot and violin graphs were created in R 3.3.4 (R Core Team 2015) with package ggplot2 (Wickham 2009) and 
show the differences between the width and height of ascospores of particular species and their hybrids. Different letters above the plot indicate significant 
difference (P < 0.05) in the size of the ascospores based on Tukey's HSD test. Boxplots show median, interquartile range, values within + 1.5 of interquartile 
range (whiskers) and outliers. 
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x ` A 2 
Fig.9 Ascospore morphology of interspecific hybrids between A. udagawae and other species. a—g. Hybrid of A. udagawae CMF ISB 2190 x A. acrensis 
IFM 57290; a-b. ascospores in light microscopy; c—g. ascospores in scanning electron microscopy; h-r. hybrid of A. udagawae CCF 4479 x A. felis NRRL 
62901; h—k. ascospores in light microscopy; I—r. ascospores in scanning electron microscopy; s—v. hybrid of A. udagawae IFM 469721 x A. wyomingensis 
CCF 4411; s—t. ascospores in light microscopy; u—v. ascospores in scanning electron microscopy. — Scale bars: a—b, h—k, s—t = 5 um; c—g, l-r, u-v = 2 um. 


162 


Persoonia — Volume 41, 2018 


CCF 5617 


CCF 5627@ 
CCF 5625@ 


IFM 60053 


A. felis 


CCF 4148 
NRRL 62901 


CCF 5610 


NRRL 62904" 


A. pseudoviridinutans IFM 59502 


CCF 5632 
CCF 4169 


A. wyomingensis [ 
CCF 44177 


' IFM 572917 
A. acrensis 


CCF 4961 
IFM 61334" 
A. arcoverdensis IFM 59922 
CBS 121595 


A. viridinutans IFM 470457 ê 


A. frankstonensis 


O 


Ascospore body width (um) 
+ oa 
d <T Q) 
| < IL r 5 


O 
' 


m A. wyomingensis 


Ascospore body height (um) 


a b 
5- 
4- 


[J A. felis 


e@elFM 57289 
eCBS 458.75 


© CCF 4416 


@IFM 57290 


NRRL 62903 


CCF 5622 
eCCF 5609 
CCF 4498 
NRRL 62900 


IMI 182127 


CCF 4419 
IMI 133982 
CCF 4959 
CCF 4960 


IFM 61345 
IFM 61338 
IFM 61362 


CBS 142233" 
CBS 142234 


! 


BJ A. wyomingensis x A. felis ë A. felis x A. pseudoviridinutans 


Fig. 10 Schematic depiction of results of interspecific mating experiments between opposite mating type isolates of heterothallic members of Aspergillus 
viridinutans species complex except of A. udagawae. Only successful mating experiments are displayed by coloured connecting lines between opposite mat- 
ing type isolates (different colours correspond to hybrids between different species); grey dashed lines indicate production of infertile ascomata; remaining 
mating experiments were negative. Boxplot and violin graphs were created in R 3.3.4 (R Core Team 2015) with package ggplot2 (Wickham 2009) and show 
the differences between the width and height of ascospores of particular species and their hybrids. Different letters above the plot indicate significant differ- 
ence (P < 0.05) in the size of the ascospores based on Tukey’s HSD test. Boxplots show median, interquartile range, values within + 1.5 of interquartile range 


(whiskers) and outliers. 
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Fig. 11 Ascospore morphology of interspecific hybrids between A. felis, A. pseudoviridinutans and A. wyomingensis. a—e. Hybrid of A. felis x A. pseudovir- 
idinutans; a—c. ascospores of hybrid CCF 4497 x IFM 59502 in light microscopy; d—e. ascospores in scanning electron microscopy: CCF 4497 x IFM 59502 
(d), CCF 4171 x IFM 59502 (e); f—I. hybrid of A. felis NRRL 62900 x A. wyomingensis CCF 4169; f—g. ascospores in light microscopy; h—I. ascospores in 


scanning electron microscopy. — Scale bars: a—c, f—g = 5 um; d—e, h-I = 2 um. 


TAXONOMY 


Aspergillus acrensis Hubka, A. Novakova, Yaguchi, Matsuz. 
& Y. Horie, sp. nov. — MycoBank MB822542; Fig. 14 


Etymology. Named after the region of origin of the ex-type strain — state 
Acre located in the northern Brazil. 


Mycelium composed of hyaline, branched, septate, smooth- 
walled hyphae. Conidial heads greyish green, loosely columnar, 
up to 140 um long, 15—25 um diam. Conidiophores uniseriate, 
arising from aerial hyphae or the basal mycelium, hyaline to 
pale yellowish brown, frequently nodding, smooth, 150—600 um 
long; stipes 3—5.5(—8) um wide in the middle; vesicles hyaline to 
greyish green, pyriform, subclavate to clavate, (6—)9—16(—20) 
um diam; phialides ampulliform, hyaline to greyish green, 4.5— 
6(-7.5) x 1.5—2.5(—3) um, covering approximately the apical 
half of the vesicle. Conidia hyaline to greyish green, globose, 
subglobose to broadly ellipsoidal, smooth-walled to delicately 
roughened, microtuberculate in SEM, 2.5-3 x 2—2.5 um (mean 
+ standard deviation, 2.8 + 0.2 x 2.4 + 0.2; length/width ratio 
1.1-1.3, 1.2 + 0.1). Heterothallic, sexual morph unknown. 


Culture characteristics (7 d at 25 °C, unless otherwise stated) 
— Colonies on MEA attained 51—62 mm diam, sparsely lanose, 
slightly raised, flat, yellowish white (ISCC—NBS No. 92) to 
pale green (No. 149), no exudate, soluble pigment light grey- 
ish yellow (No. 101), reverse light greenish yellow (No. 101) 
to brilliant greenish yellow (No. 98). Colonies on CYA attained 
33-48 mm diam, floccose, slightly raised, flat to slightly radi- 
ally furrowed, yellowish white (No. 92) to greenish white (No. 
153), sporulation in the colony centre pale green (No. 149) to 
greyish green (No. 150), no exudate, soluble pigment dark 
greyish yellow (No. 91), reverse deep yellow (No. 85), light 
olive brown (No. 94) to moderate olive brown (No. 95) with 
light yellow (No. 86) margin. Colonies on CYA at 37 °C grow 
more rapidly compared to 25 °C and attained 60—70 mm diam, 
lanose, slightly raised, flat to radially furrowed, white mycelium 
in margins, sporulation light olive grey (No. 112) to olive grey 
(No. 113), no exudate, no soluble pigment, reverse colourless, 
moderate yellow (No. 87) to greyish yellow (No. 90). Colonies 
on CZA attained 36—48 mm diam, lanose, slightly raised, flat, 
yellowish white (No. 92), no exudate, no or light greyish yel- 
low (No. 101) soluble pigment, reverse light yellow (No. 86), 


164 


Persoonia — Volume 41, 2018 


ur 


Fig. 12 Sexual morph morphology of homothallic species from Aspergillus viridinutans complex. a—d. Aspergillus aureolus isolates IFM 47021" (a—b, d) and 
IFM 46584 (c); a. Macromorphology of ascomata after 3 wk of incubation on MEA at 37 °C; b. ascospores in light microscopy; c—d. ascospores in scanning 
electron microscopy; e—h. Aspergillus siamensis isolate IFM 59793"; e. macromorphology of ascomata after 3 wk of incubation on MEA at 37 °C; f. ascospores 
in light microscopy; g—h. ascospores in scanning electron microscopy. — Scale bars: b, f = 5 um; c—d, g—h = 2 um. 


light greenish yellow (101) to brilliant greenish yellow (No. 98). 
Colonies on YES lanose, yellowish white (No. 92), irregularly 
furrowed, no exudate, soluble pigment brilliant yellow (No. 83), 
reverse brilliant yellow (No. 83). Colonies on CY20S attained 
58-65 mm diam, lanose, slightly raised, flat, yellowish white 
(No. 92), no exudate, no soluble pigment, reverse moderate 
brown (No. 58) to moderate reddish brown (No. 43). Colonies 
on CREA attained 32-35 mm diam, sparsely lanose, plane, 
mycelium yellowish white, no visible sporulation, reverse strong 
brown (No. 55), no acid production. Growth on MEA at 45 °C, 
no growth on MEA at 47 °C. 

Exometabolites — Isolate IFM 57291 produced an aszona- 
pyrone, a fumigatin, tryptoquivalines, tryptoquivalones; isolate 
IFM 57290 an aszonapyrone, fumagillin, fumigatins, helvolic 
acid, pseurotin A, tryptoquivalines, and a trytoquivalone; iso- 
late CCF 4959 pseurotin A, viriditoxin and several potential 
naphtho-gamma-pyrones; CCF 4960 antafumicins, fumagillin, 
a fumigatin, helvolic acid, pseurotin A, and a tryptoquivalone; 
and CCF 4961 an aszonapyrone, fumagillin, fumigatins, 
pseurotin A, tryptoquivalines and tryptoquivalones. In general, 
similar metabolites are also produced by the two most closely 
related species, i.e., A. aureolus and A. udagawae. Aspergillus 
aureolus produces fumagillin, helvolic acid, pseurotin A, trypto- 


quivalines, tryptoquivalones and viriditoxin as well as several 
unique yellow secondary metabolites. Aspergillus udagawae 
produces fumagillin, fumigatins, tryptoquivalines and trypto- 
quivalones (Frisvad & Larsen 2015a). 


Specimens examined. Brazit, State of Acre, Xapuri, grassland soil in 
a cattle farm, 6 Nov. 2001, Y. Horie (holotype IFM 57291H, isotypes PRM 
935088 and PRM 935089, culture ex-type IFM 57291" = CCF 4670'); State of 
Amazonas, Manaus, forest soil in tropical rain forest, 11 Nov. 2001, Y. Horie, 
culture IFM 57290 (= CCF 4666). — Romania, Movile cave, above the Lake 
Room, cave sediment, 8 June 2014, A. Novakova, culture CCF 4959; Movile 
cave, cave corridor, cave sediment, 8 June 2014, A. Novakova, culture CCF 
4960; Movile cave, Lake Room, cave sediment, 8 June 2014, A. Novakova, 
culture CCF 4961. 


Notes — The morphology of A. acrensis strongly resembles 
that of several other A. viridinutans complex members. The 
closely related taxa A. aureolus and A. siamensis are readily 
distinguished from A. acrensis by the production of ascomata 
under standard cultivation conditions (both are homothallic). 
Aspergillus viridinutans and A. frankstonensis grow more slowly 
at 25 °C and have smaller vesicles. The macromorphology of 
colonies and micromorphology of the asexual morph does not 
distinguish A. acrensis reliably from A. arcoverdensis, A. felis, 
A. pseudoviridinutans, A. udagawae and A. wyomingensis. 
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Fig. 13 Conidia with micro-tuberculate surface ornamentation pattern observed by scanning electron microscopy. a. Aspergillus acrensis IFM 57290; 
b. A. arcoverdensis IFM 613347; c. A. aureolus IFM 46584; d. A. felis CBS 1302457; e. A. felis NRRL 62900 (ex-type of A. parafelis); f. A. felis NRRL 62903 
(ex-type of A. pseudofelis); g. A. frankstonensis CBS 142234; h. A. pseudoviridinutans CBS 458.75; i. A. siamensis IFM 59793"; j. A. udagawae IFM 469727; 
k. A. viridinutans IFM 47045"; |. A. wyomingensis CCF 4414. — Scale bars = 2 um. 
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Fig. 14 Micromorphology and macromorphology of Aspergillus acrensis. a—e. Colonies of IFM 57291" incubated 7 d at 25 °C on MEA, CYA, CZA, YES, and 
on CYA at 37 °C (from left to right); f—j. reverse of colonies of IFM 57291" incubated 7 d at 25 °C on MEA, CYA, CZA, YES, and on CYA at 37 °C (from left to 
right); k—n. conidiophores; o. conidia. — Scale bars = 10 um. 
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Some of these species can be differentiated each from the 
other by their characteristic sexual morph, but the production 
of ascomata was not induced in A. acrensis despite our at- 
tempts, similarly to A. arcoverdensis and A. pseudoviridinutans. 
Although isolate IFM 57290 was successfully crossed with 
isolates of A. udagawae IFM 46972 and CMF ISB 2190 in vitro, 
both the width and height of ascospores were statistically dif- 
ferent from A. udagawae. Also, abnormalities in the shape and 
superficial ornamentation (Fig. 9) were present in a significant 
number of spores (equatorial crests were absent in ~ 50 % of 
ascospores). Reliable identification of A. acrensis can currently 
only be achieved by molecular methods. 


Aspergillus udagawae Horie et al., Mycoscience 36: 199. 
1995. 


Epitypification. Brazi., São Paulo State, Botucatu, Lagoa Seka Avea, soil 
in a plantation, 23 Aug. 1993, M. Takada (holotype CBM-FA-0711, designated 
by Horie et al. (1995), epitype designated here PRM 945579, isoepitypes 
PRM 945580 and 945581, MycoBank MBT378451, culture ex-epitype IFM 
46972 = CBS 114217 = DTO 157-D7 = CBM-FA 0702 = KACC 41155 = CCF 
4558). 


Notes — Horie et al. (1995) designated the specimen CBM- 
FA-0711 as a holotype of A. udagawae, a dried culture with 
ascomata created by crossing the isolates CBM-FA-0702 
(MAT1-1-1) x CBM-FA-0703 (MAT1-2-1). Although this speci- 
men demonstrates the sexual and asexual morph of the life 
cycle, it is not suitable for the purposes of the recent taxonomy 
for several reasons. First of all, it is not clear which of the two 
cultures contained within the type should be considered the 
ex-holotype culture. Additionally, interspecific hybrids can be 
induced by crossing opposite mating type strains of unrelated 
species in vitro as shown in this study and some previous 
studies (see Discussion), and deposition of a resultant ‘hybrid’ 
type could lead to ambiguities. Although this second argument 
does not apply to A. udagawae as both isolates included in the 
holotype are closely related phylogenetically, we believe that 
a more clearly defined type of this species will facilitate future 
taxonomic work. Because it is not possible to recognize which 
part of the holotype belongs to particular isolate, lectotype 
designation (in this case part of holotype specimen) is difficult. 
For this reason we decided to select an epitype PRM 945579 
derived from the IFM 46972 (= CBM-FA 0702) culture. 


DISCUSSION 


Changing species concepts in the AVSC 


The AVSC members show considerable phenotypic variability 
but usually share production of nodding heads (some vesicles 
borne at an angle to the stipe) and relatively poor sporulation 
with abundant aerial mycelium. All species have a maximum 
growth temperature of 42 or 45 °C and the macromorphology 


and diameter of their colonies are similar, except for A. viridi- 
nutans and A. frankstonensis, which grow more slowly than 
remaining species. In addition, the morphology of conidiophores 
and conidia is relatively uniform across species, including the 
superficial ornamentation of conidia as shown here (Fig. 13). 
For these reasons heterothallic AVSC members have resisted 
taxonomic classification and were only identified to a species 
complex level, until recently. 


Due to the absence of taxonomically informative characters, 
most recently described species in the AVSC were delimited 
using the GCPSR rules. Using this approach, the species are 
recognized based on concordance between single-gene phy- 
logenies and the absence of tree incongruities. The GCPSR 
has found wide application in the taxonomy of fungi (Dettman 
et al. 2006, Hubka et al. 2013a, Peterson et al. 2015, Visagie 
et al. 2017). Huge progress has been made recently in the 
development of statistical methods for multilocus species deli- 
mitation, driven by advances in the multispecies coalescent 
model (Bouckaert et al. 2014, Flot 2015, Fontaneto et al. 
2015, Schwarzfeld & Sperling 2015, Jones 2017). Although 
the ideology of MSC delimitation methods is relatively similar to 
GCPSR, these methods are more robust because the species 
are delimited in three steps, i.e., species discovery, species 
tree construction and species validation (Carstens et al. 2013). 
The determination of species boundaries is more objective in 
contrast to GCPSR rules that are based on relatively subjective 
evaluation and comparison of single-gene trees. In addition, 
MSC methods are able to deal better with phenomena such 
as incomplete lineage sorting, recombination or non-reciprocal 
monophyly that lead to incongruences between single-gene 
trees. Compared to the phylogenetic analysis of concatenated 
gene datasets (including partitioned datasets) and in part also 
GCPSR, the MSC methods are less prone to over-delimitation 
of species (Degnan & Rosenberg 2006, Kubatko & Degnan 
2007, Heled & Drummond 2010, Rosenberg 2013), especially 
when the results of multiple delimitation methods are compared 
in one analysis. 


The GCPSR rules together with evaluation of limited phenotypic 
data were recently used for description of A. felis, A. arcover- 
densis and A. frankstonensis in the AVSC (Barrs et al. 2013, 
Matsuzawa et al. 2015, Talbot et al. 2017). Genealogical analy- 
sis using five genetic loci was carried out for delimitation of 
A. parafelis, A. pseudofelis and A. pseudoviridinutans, three 
close relatives of A. felis (Sugui et al. 2014). Although the 
authors found no conflict between single-gene phylogenies, 
only two isolates of each of these four species were used in 
analysis, and sequences of A. felis, A. parafelis and A. pseudo- 
felis strains included were almost invariable. These isolates 
did not cover sufficiently the genetic diversity of these species 
as shown here. Species delimitation results based on MSC 
in this study showed that A. parafelis and A. pseudofelis are 
included in the genetically diverse lineage of A. felis (Fig. 3). 


Table 6 Highest intraspecific pairwise genetic distances in members of Aspergillus viridinutans complex (%). 


Highest genetic distances between two isolates according to different genetic loci 


Species (no. of isolates) 


benA CaM RPB2 act mcm7 tsr1 
A. acrensis (5) 0.2 0.9 0.2 0 ND ND 
A. arcoverdensis (13) 0 0.9 0.5 1.4 ND ND 
A. aureolus (4) 0.4 0 0.1 0 ND ND 
A. felis (35) 4.2 2.4 0.6 25 1.3 3.3 
A. frankstonensis (2) 0 0.2 0 0 0 ND 
A. pseudoviridinutans (8) 2.6 2.2 1.9 2.1 0.7 1.4 
A. siamensis (2) 0 0.1 0.1 0 ND ND 
A. udagawae (25) 1.1 2.8 1.2 4.9 ND ND 
A. wyomingensis (15) 0.4 0.9 0.4 0.9 ND ND 


ND, not determined. 
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The intraspecific pairwise genetic distances in A. felis (Table 6) 
range from 0.6 % (RPB2) to 4.2 % (benA). Similarly, pairwise 
genetic distances in A. udagawae (Table 6) are 1.1 % (benA) 
to 4.9 % (act). Such high intraspecific diversity in these genetic 
loci is unusual in Aspergillus and it reflects the intense recom- 
bination. Thus, when only limited number of strains from such 
species are selected for phylogenetic analysis, the results of 
species delimitation techniques may be biased and prone to 
overestimate the number of species. As we have shown here, 
this was clearly the case in the study of Sugui et al. (2014). 
This problem is probably widespread in current fungal taxonomy 
and limits possibilities of correct species boundaries delimita- 
tion. Also in this study, the number of strains of some closely 
related and phenotypically similar species is underrepresented, 
e.g., A. viridinutans and A. frankstonensis. In these cases, the 
species boundaries cannot be reliably defined using neither 
GCPSR rules nor MSC-based methods used in this study. 


Clinically relevant species and their identification in clini- 
cal setting 


Although sect. Fumigati harbours many important pathogenic 
species, members of the AVSC have been overlooked by both 
clinicians and mycologists until recently. The presence of these 
soil-borne species in clinical material was first reported by Katz 
et al. (2005) who examined phylogenetic positions of several 
‘atypical’ (poorly sporulating) clinical isolates of A. fumigatus. 
The majority of these strains grouped with, but were not identi- 
cal to, A. viridinutans and A. aureolus from the AVSC. Since 
then many similar epidemiological and clinical studies have 
reported the pathogenic role of AVSC members in humans and 
animals, as reviewed by Talbot & Barrs (2018). In humans the 
most common manifestation of disease is chronic invasive pul- 
monary aspergillosis in immunocompromised patients. AVSC 
species are also frequently reported as a cause of sino-orbital 
aspergillosis (SOA) in cats that is chronic, but frequently fatal. 
In contrast to humans and dogs, the disease usually affects 
ostensibly immunocompetent cats. This increasingly recognised 
clinical entity is most frequently caused by the AVSC species 
and less frequently by other cryptic species in sect. Fumigati 
(Barrs et al. 2012, 2013, 2014). 


Based on the species boundaries redefined in this study, the 
AVSC encompasses four species that are confirmed oppor- 
tunistic pathogens. According to a number of reported cases, 
in humans A. udagawae is the most important opportunistic 
pathogenic from the AVSC followed by A. felis and A. pseudo- 
viridinutans. In contrast, SOA in cats is most commonly caused 
by A. felis and much less frequently by A. udagawae and 
A. wyomingensis (Barrs et al. 2013, 2014). 


Medically important species from the AVSC demonstrate ele- 
vated minimum inhibitory concentrations (MICs) of itraconazole 
and voriconazole in vitro, and a variable susceptibility to ampho- 
tericin B, while posaconazole and echinocandins have potent 
in vitro activities (Lyskova et al. 2018). Since the intraspecific 
variation in MICs of particular antifungals is usually high, the use 
of reliable methods for MIC determinations takes precedence 
over correct identification to a species level. The latter may be 
challenging or even impossible in the clinical setting. However, 
identification to the level of species complex and differentia- 
tion from A. fumigatus is important due to strikingly different 
antifungal susceptibility patterns. 


In contrast to A. fumigatus, the AVSC species do not grow 
at 47 and 50 °C, usually sporulate less and a proportion of 
their vesicles are borne at an angle to the stipe. In addition, 
some isolates produce acidic compounds detectable on CREA 
(Barrs et al. 2013, Novakova et al. 2014, Talbot & Barrs 2018). 
Despite the fact that ITS rDNA region sequences are not avail- 
able for all AVSC members, this universal marker for fungal 


species identification and barcoding can be used to achieve 
identification to a species complex level. The sequences from 
all six protein-coding genes included in this study (Table 2) have 
sufficient discriminatory power for species level identification of 
all clinically relevant species. Among these genes, sequences of 
B-tubulin and calmodulin belong to the most commonly used in 
the clinical practice when correct identification is required (epi- 
demiological studies, outbreak investigations or when dealing 
with infections refractory to antifungal therapy). However, the 
discrimination between A. felis and A. pseudoviridinutans can 
be limited when using the B-tubulin gene due to the incomplete 
lineage sorting phenomenon detected in this study (Fig. 3). 


Additionally, the increasingly used method of matrix-assisted 
laser desorption/ionization time-of-flight mass spectrometry 
(MALDI-TOF MS) gives promising results for rapid and accurate 
discrimination between A. fumigatus and other clinically relevant 
aspergilli from sect. Fumigati (Alanio et al. 2011, Nakamura et 
al. 2017). The development of more robust, curated and ac- 
cessible MALDI-TOF spectrum databases should enable the 
implementation of MALDI-TOF MS for routine identification of 
less common aspergilli in future. Several PCR assays targeting 
protein-coding or microsatellite loci have also been developed 
and show good efficiency in discrimination of less common path- 
ogenic species in sect. Fumigati (Araujo et al. 2012, Fernandez- 
Molina et al. 2014, Chong et al. 2017). 


Mating behaviour in the AVSC — heterothallic species 


The increasing availability of PCR-based tools for identification 
of fungal genes responsible for sexual and somatic incompati- 
bility has facilitated the ability to induce the sexual morph in 
fungi (Dyer & O’Gorman 2011). The characterisation of mating 
type (MAT) genes became routine when inducing the sexual 
morph of heterothallic species in vitro. Using this approach, 
the sexual morph has been induced recently in at least five 
members of sect. Fumigati (O’Gorman et al. 2009, Barrs et 
al. 2013, Swilaiman et al. 2013, Novakova et al. 2014, Hubka 
et al. 2017). The discovery of a sexual cycle in pathogenic 
and mycotoxigenic fungi has many important consequences, 
because fungi with a functional sexual cycle have greater po- 
tential to increase their virulence and to develop resistance to 
antifungals, fungicides, etc. (Kwon-Chung & Sugui 2009, Lee 
et al. 2010, Swilaiman et al. 2013). 


Here, we induced the sexual morph with a relatively high rate of 
success in A. felis, A. udagawae and A. wyomingensis (Fig. 6). 
We demonstrated that ascospores of these three species have 
relatively stable morphology (Fig. 7) and that the size of their 
ascospores is significantly different from one another (Fig. 6) 
and can be differentiated by equatorial crest length (Table 5). 
However, not all opposite mating type strains of the same spe- 
cies are able to produce ascomata in vitro as demonstrated in 
all three mentioned species (Fig. 6). A similar decline in mating 
capacity was also demonstrated in previous studies on the 
AVSC (Sugui et al. 2010, Novakova et al. 2014), but also in 
A. lentulus (Swilaiman et al. 2013) and A. fumigatus (O’Gorman 
et al. 2009). These species require relatively rigid conditions 
to complete their sexual cycle and some crosses produce low 
numbers of or infertile ascomata or do not mate at all (Balajee 
et al. 2006, Yaguchi et al. 2007, Kwon-Chung & Sugui 2009, 
Sugui et al. 2010, Novakova et al. 2014). For instance, fertility 
between two opposite mating-type isolates may be influenced 
by the vegetative incompatibility genes (Olarte et al. 2015), 
regulators of cleistothecium development and hyphal fusion 
(Szewczyk & Krappmann 2010). 


We were not able to induce the sexual morph in three heterothal- 
lic members of the AVSC, i.e., A. acrensis, A. arcoverdensis 
and A. pseudoviridinutans, despite the relatively high number of 
opposite mating-type strains that was available for the mating 
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Table 7 Genetic similarities between the ex-type isolates of Aspergillus viridinutans complex members based on identities from BLAST similarity search’. 


Genetic similarities between species: benA/ CaM I RPB2 (%) 


Species 


10. 


A. acrensis 


1. 


94.5/95.2/98.0 


2. A. arcoverdensis 


94.5/95.6/98.1 


99.6/98.8/99.0 


3. A. aureolus 


93.4/96.8/97.6 92.4/95.9/97.8 


92.0/95.6/97.7 


4. A. felis 


95.6/97.1/98.3 95.3/94.9/98.2 92.6/96.2/97.7 


95.3/94.7/98.0 


5. A. frankstonensis 


94.7/95.2/97.6 95.7/96.0/97.4 94.9/95.5/97.8 95.5/97.6/98.1 96.0/95.3/97.5 


6. A. pseudoviridinutans 


95.5/95.6/98.5 96.7/96.0/98.9 93.0/95.7/98.2 95.6/94.7/98.0 95.5/95.4/97.9 


96.6/95.8/98.9 


7. A. siamensis 


94.7/95.6/98.2 97.4/97.1/99.1 92.0/95.9/97.9 95.3/95.1/98.1 94.5/95.6/97.7 96.2/96.3/99.1 


97.4/96.8/99.0 


8. A. udagawae 


96.5/97.3/99.1 95.5/95.1/98.6 93.8/95.4/98.2 97.5/97.8/98.8 96.5/94.7/97.9 96.3/95.3/98.8 95.6/95.3/98.8 


95.3/94.8/98.6 


9. A. viridinutans 


94.5/96.0/97.8 96.0/96.5/98.3 92.1/96.3/97.5 95.4/95.8/97.6 94.9/95.9/97.3 96.9/96.5/98.9 95.8/96.9/98.6 95.6/95.7/98.3 


95.8/96.5/98.6 


1 nucleotide BLAST with default setting (http://blast.ncbi.nlm.nih.gov/Blast.cgi). 


10. A. wyomingensis 


assays (Fig. 6). It is not clear if these species require different 
conditions for successful mating, if there are other unidentified 
pre-zygotic mating barriers between opposite mating type 
strains, or if they have lost the ability to complete their sexual 
cycle. The evidence that two of these species were able to mate 
with different species from AVSC makes the last possibility im- 
probable (Fig. 8, 10). These hybrids can be differentiated from 
A. udagawae and A. felis, respectively, by their dimensions (Fig. 
8, 10) and surface ornamentation (Fig. 9, 11; Table 5). It demon- 
strates that both A. acrensis and A. pseudoviridinutans should 
be treated as separate taxonomic entities from their related 
species. Similar deviations in size and surface ornamentation 
of ascospores were demonstrated in other interspecific hybrids 
(Fig. 8—11) when they were compared to parental species. 


Mating behaviour in the AVSC — homothallic species 


Although homothallic species prevail over heterothallic in sect. 
Fumigati (Fig. 1), only two homothallic species are present 
in the AVSC. It is supposed that heterothallism is ancestral 
to homothallism in fungi (Nauta & Hoekstra 1992), including 
Aspergillus (Rydholm et al. 2007, Lee et al. 2010). It is obvious 
from phylogenetic studies across different subgenera of Asper- 
gillus, that reproductive strategy is evolutionary conservative 
and homothallic as well as heterothallic (or asexual) species are 
typically clustered in clades with a uniform reproductive strategy. 
For instance in subg. Aspergillus, the 31 currently accepted 
species of sect. Aspergillus are all homothallic (Chen et al. 
2016a) while sister sect. Restrictiencompasses 20 asexual and 
only one distantly related homothallic species, A. halophilicus 
(Sklenář et al. 2017). Similarly, subg. Polypaecilum harbours 
only asexual species (Martinelli et al. 2017, Tanney et al. 2017). 
Asexual species also predominate in subg. Circumdati (Jurjević 
et al. 2015) although most, if not all, probably have a cryptic 
sexual cycle as highlighted by sexual morph induction in A. 
flavus, A. nomius, A. parasiticus, A. terreus and A. tubingen- 
sis (Horn et al. 2009a, b, 2011, 2013, Arabatzis & Velegraki 
2013). A strikingly different situation is present in subgenera 
Nidulantes (Chen et al. 2016a, Hubka et al. 2016a), Fumigati 
(Fig. 1) and Creme (Hubka et al. 2016b) where heterothallic 
and homothallic species interchange like a mosaic along the 
phylogenetic tree. 


Common genetic distances between closely related sister spe- 
cies across aspergilli usually range between 2—4 % in benA 
and CaM loci and 1—2 % in RPB2 locus; the situation in AVSC is 
very similar (Table 7). Interestingly, there are only few examples 
of closely related homothallic and heterothallic/asexual species 
in Aspergillus despite their common origin. Genetic similari- 
ties between related couples of homothallic and heterothallic/ 
asexual exceeding 95 % are rare, with only two examples in 
subg. Circumdati and one in subg. Cremei (Table 8). Section 
Fumigati is exceptional because it contains at least five pairs 
of highly related homothallic and heterothallic species (Table 8; 
Fig. 1). Aspergillus acrensis, described here, and A. aureolus 
represent the most closely related pair across genus Aspergil- 
lus (Table 8) and thus could be an ideal model for studying the 
evolution of reproductive modes. If we accept the hypothesis 
about the derived origin of homothallic species, it is probable 
that A. aureolus evolved in the lineage of A. acrensis rela- 
tively recently, due to the extremely low genetic distances of 
both species. This is also likely the reason why the multilocus 
species delimitation method STACEY and also some single- 
locus methods failed to segregate A. acrensis from A. aureolus 
(Fig. 2) in this study. 


Interspecific hybridization in fungi and its consequences 


Interspecific hybridization is an important process affecting spe- 
ciation and adaptation of micro- and macroorganisms, however, 
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Table 8 Genetic similarities between selected homothallic species and their most closely related heterothallic / anamorphic relatives across diversity of the 


genus Aspergillus. 


Homothallic species (section) — closest heterothallic / anamorphic species 


Genetic similarities (%): 


benA/ CaM / RPB2' 
subg. Aspergillus 
A. halophilicus (Restricti) — any species < 89 
A. montevidensis (Aspergillus) — any species < 88 


subg. Circumdati 
A. alliaceus (Flavi) — A. lanosus 
A. muricatus (Circumdati) — A. ochraceus 
A. neoflavipes (Flavipedes) — A. micronesiensis 
A. neoniveus (Terrei) — any species 


subg. Cremei 
A. chrysellus (Cremei) — A. wentii 
A. cremeus (Cremei) — any species 
A. stromatoides (Cremei) — any species 


subg. Fumigati 

acanthosporus (Clavati) — A. clavatus 
aureolus (Fumigati) — A. acrensis 

cejpii (Clavati) — any species 

fischeri (Fumigati) — A. fumigatus 
posadasensis (Clavati) — A. clavatus 
quadricinctus (Fumigati) — A. duricaulis 
siamensis (Fumigati) — A. wyomingensis 
waksmanii (Fumigati) — A. nishimurae 


subg. Nidulantes 

falconensis (Nidulantes, A. nidulans clade) — A. recurvatus 
monodii (Usti) — any species 

nidulans (Nidulantes, A. nidulans clade) — any species 


purpureus (Nidulantes, A. spelunceus clade) — any species 
undulatus (Nidulantes, A. stellatus clade) — any species 


DS>>BS>EE >>> >>> 


discophorus (Nidulantes, A. aeneus clade) — A. karnatakaensis 


pluriseminatus (Nidulantes, A. multicolor clade) — any species 


96.4 / 95.7 / 99.1 
<91 
94.8 / 91.9 / 97.5 
< 90 


97.1/97.2/97.7 
<91 
<93 


<93 

99.6 / 98.8 / 99.0 
< 88 

94.3 / 94.5 / 97.9 
95.1 / 92.6 / 93.5 
92.6 / 95.0 / 99.1 
97.1 / 96.5 / 98.9 
97.8 / 98.4 / 96.6 


< 92 
£93 
< 90 
< 92 
< 92 
< 90 
< 89 


1 If none of three genetic similarities exceed 95 %, the values are replaced by only one highest value (usually RPB2 locus). 


relatively little is still Known about the frequency of hybridiza- 
tion in fungi and its role in evolution of fungal species. Fungal 
hybrids may form either by a partial or complete sexual cycle 
or by a parasexual process. Mating between two species may 
be prevented by pre-zygotic barriers (e.g., gamete recognition) 
and various post-zygotic barriers (developmental problems, 
hybrid viability and ability to reproduce, etc.). The disagreement 
between phylogenetic/morphological species concepts and bio- 
logical species compatibilities has been repeatedly described 
in fungi. Phylogenetic divergence in some fungal groups might 
have preceded development of reproductive barriers as shown 
by interspecific hybrids induced in vitro between primary hu- 
man and animal pathogenic Trichophyton species (Kawasaki 
et al. 2009, 2010, Anzawa et al. 2010, Kawasaki 2011), op- 
portunistic pathogenic Candida albicans and C. dubliniensis 
(Pujol et al. 2004), members of Aspergillus sect. Fumigati 
(Sugui et al. 2014, Talbot et al. 2017), mycotoxigenic A. flavus 
and A. parasiticus (Olarte et al. 2015), A. flavus and A. mini- 
sclerotigenes (Damann & DeRobertis 2013), phytopathogenic 
species from the Fusarium graminearum complex (Bowden & 
Leslie 1999) and species of Neurospora (Dettman et al. 2003). 
Natural interspecific hybrids resulting from recombination be- 
tween species or parasexual reproduction are most commonly 
reported and have been extensively studied in saprophytic 
yeasts (González et al. 2008, Sipiczki 2008, Nakao et al. 2009, 
Louis et al. 2012), the plant endophyte Epichloë (Cox et al. 
2014, Charlton et al. 2014, Shymanovich et al. 2017) and in 
various plant pathogenic fungi including species of Fusarium 
graminearum complex (O'Donnell et al. 2004, Starkey et al. 
2007), Ophiostoma (Brasier et al. 1998, Solla et al. 2008), 
Microbotryum (Gladieux et al. 2010), Melampsora (Spiers & 
Hopcroft 1994, Newcombe et al. 2000), Botrytis (Staats et al. 
2005), Verticillium (Inderbitzin et al. 2011) and Heterobasidion 
(Gonthier et al. 2007, Lockman et al. 2014). 


Considering that in vitro induction of hybrids is relatively suc- 
cessful, it is surprising that reports on the isolation of naturally 
occurring hybrids are infrequent in human and animal patho- 
genic fungi. It suggests that post-zygotic mating barriers play 
a fundamental role in the maintenance of species boundaries. 
Naturally occurring hybrids have been detected in yeasts and 
dimorphic fungi, including between Candida spp. (Schröder et 
al. 2016), Malassezia spp. (Wu et al. 2015), Cryptococcus neo- 
formans and C. gattii (Bovers et al. 2006, 2008, Kwon-Chung & 
Varma 2006, Aminnejad et al. 2012) and Coccidioides immitis 
and C. posadasii (Johnson et al. 2015). However, to date, 
reports on these hybrids in filamentous fungi are restricted to 
the Neocosmospora solani complex (Short et al. 2013, 2014). 
Species definition has become a controversial issue in some 
of these species complexes with naturally occurring hybrids 
because of differing opinions on species concepts among 
taxonomists (Kwon-Chung & Varma 2006, Kawasaki 2011, 
Kwon-Chung et al. 2017). 


Even in cases where interspecific hybrids with high fitness and 
fertility can be demonstrated, the intensity of gene flow between 
natural populations must be sufficient to oppose genetic drift in 
order to have a significant impact on genetic isolation of species. 
In fungi, these processes cannot be evaluated rigorously by in 
vitro mating assays, as these cannot be extrapolated fully to a 
natural setting (Starkey et al. 2007, Sugui et al. 2014, Hubka et 
al. 2015a). Indeed, natural interspecific hybrids have never been 
reported for the majority of species that readily produce hybrids 
in vitro, including Aspergillus, dermatophytes and Neurospora. 
The MSC and GCPSR approaches provide practical tools for 
evaluating the significance of gene flow between natural popu- 
lations and for assessing species limits. The interpretation 
of in vitro mating assays without a robust phylogeny is thus 
controversial, because a number of clearly phylogenetically, 
morphologically and ecologically distinct species lack effective 
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reproductive barriers. In addition, the evaluation of biological 
species limits using mating assays requires determination of 
the fitness and fertility of progeny, which is demanding in both 
time and cost. 


In general, mating success between different species under 
laboratory conditions is much lower compared to intraspecific 
mating, suggesting strong reproductive isolation between 
species and adherence to the biological species concept. In 
agreement with this, only a limited number of strains with ex- 
ceptional mating capacity are usually capable of interspecific 
hybridization with strains of different species, e.g., A. udagawae 
strain IFM 46972 (Fig. 8) or A. pseudoviridinutans strain IFM 
59502 (Fig. 10). 


Several studies demonstrated that interspecific hybrids express 
genetic abnormalities or have decreased fertility and viability. 
Genetic analysis of the progeny of a cross between F. asiaticum 
x F. graminearum detected multiple abnormalities that were 
absent in intraspecific crosses of F. graminearum, i.e., pro- 
nounced segregation distortion, chromosomal inversions, and 
recombination in several studied linkage groups (Jurgenson et 
al. 2002, Gale et al. 2005). Matings between C. neoformans x 
C. gattii produced only a low percentage of viable progeny. It 
has been suggested that C. neoformans and C. gattii produce 
only stable diploid hybrids, but not true recombinants (Kwon- 
Chung & Varma 2006). Although Olarte et al. (2015) obtained 
hybrid progeny of A. flavus and A. parasiticus, fertile crosses 
were rare and involved only one parental strain of A. flavus. 
Viable ascospores were extremely rare, suggesting extensive 
genetic incompatibility and post-zygotic incompatibility mecha- 
nisms. Morphologically, the progeny differed from parental 
strains in growth rate, sclerotium production, stipe length, co- 
nidial head seriation and conidial features (Olarte et al. 2015). 
Decreased viability of hybrid ascospores was also detected 
among Neurospora spp. (Dettman et al. 2003) and in Aspergil- 
lus sect. Fumigati, in addition to abnormalities in their surface 
ornamentation visualised by SEM (Sugui et al. 2014), which is 
in agreement with the present study (Fig. 9, 11). Apart from as- 
cospore ornamentation, we also found significant differences in 
hybrid ascospore dimensions from parental species (Fig. 8, 10). 


The relatively recent globalization of trade in horticultural and 
agricultural plants, and introduction of non-native plant species 
has resulted in the inadvertent introduction of alien plant patho- 
gens into non-endemic areas, contributing to the emergence of 
some devastating plant diseases (Brasier 2001, Mehrabi et al. 
2011, Dickie et al. 2017). Anthropogenic activities or changes 
in the distribution of fungi (e.g., in response to climate changes) 
may bring together related, previously allopatric pathogenic 
species. Subsequent interspecific hybridization could give rise 
to pathogens with new features, including adaptation to new 
niches and host species, and varying degrees of virulence, as 
evidenced in Verticillium longisporum, Zymoseptoria pseudo- 
tritici, Blumeria graminis f. sp. triticale, and hybrids between 
Ophiostoma novo-ulmi and O. ulmi (Brasier 2001, Schardl & 
Craven 2003, Depotter et al. 2016). 


As far as we know, the occurrence of Aspergillus interspecific 
hybrids in nature has not been proven despite successful hybri- 
dization of some species in vitro. However, there is no reason 
to assume that this phenomenon does not occur occasionally. 
Genetic recombination similar to that found in intraspecific mat- 
ing occurred in half of the progeny produced by mating A. fumi- 
gatus with A. felis, while the other half were probably diploids or 
aneuploids (Sugui et al. 2014). Progeny resulting from mating 
between A. flavus and A. minisclerotigenes was fertile when 
crossed with parental strains and the frequency of successful 
matings was similar to that within pairs of A. flavus and A. mini- 
sclerotigenes strains, respectively (Damann & DeRobertis 


2013). Ultimately, the viable hybrid must present some charac- 
teristics that promotes its survival (Turner et al. 2010, Mixao & 
Gabaldón 2018). For instance Olarte at al. (2015) showed that 
some F. progeny of A. flavus x A. parasiticus produced higher 
aflatoxin concentrations compared to midpoint parent aflatoxin 
levels, and some hybrids synthesized G aflatoxins that were not 
produced by the parents. This suggested that hybridization is 
an important diversifying force generating novel toxin profiles 
(Olarte et al. 2015). Although interspecific hybridization in as- 
pergilli is a relatively newly discovered phenomenon, it is likely 
to have played an important role in the evolution of the genus. 


The relationship between hybridization and changes in virulence 
potential is not well understood in human and animal fungal 
pathogens but its role in the emergence of novel plant fungal 
pathogens is well documented, as discussed. The evidence of 
biological compatibility between major pathogens in Aspergillus 
sect. Fumigati sheds new light on possible interspecific transfer 
of virulence genes, genes responsible for antifungal resist- 
ance, and other genes influencing adaptation of these fungi to 
a changing environment. Further studies should elucidate to 
what extent interspecific hybridization shaped the evolution of 
these opportunistic pathogens. 


CONCLUSIONS 


Based on consensus results of species delimitation methods and 
after evaluation of mating assay results and phenotypic data, 
we now recognise 10 species within the AVSC. This number 
comprises nine previously recognised and one new species 
proposed here. Aspergillus pseudofelis and A. parafelis are 
placed in synonymy with A. felis. All four genetic loci used for 
phylogenetic analysis across the AVSC have sufficient variability 
for reliable species identification and can be used as DNA bar- 
codes. Though more laborious, the MSC are a suitable tool for 
delimitation of genetically diverse cryptic species in cases where 
classical phylogenetic, morphological and mating compatibility 
data do not yield satisfactory results. 
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